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A New Approach to Space Diversity Combining 
in Microwave Digital Radio 


By Y.S. YEH and L. J. GREENSTEIN* 
(Manuscript received May 30, 1984) 


In this paper we describe a new approach to dual-channel space diversity 
combining in microwave digital radio. This approach features (1) adaptive 
control of the relative amplitudes and phases of the two branch gains; and (2) 
a search strategy, based on noncoherent spectrum measurements at the 
combiner output, that simultaneously accounts for both dispersion and noise. 
Computer programs have been developed to simulate the search process and 
to analyze the resulting performance. Eight representative channel response 
pairs are postulated and performance results are presented for each. They 
show that the scheme provides a high degree of channel equalization over 
bandwidths up to at least 40 MHz, and that, in receivers not using adaptive 
equalizers, it offers major improvements in detection performance over selec- 
tion diversity. 


I. INTRODUCTION 


In a terrestrial digital radio link, frequency selective fading caused 
by multipath propagation presents the major threat to system availa- 
bility. Efforts to reduce channel dispersion, and thus to increase 
availability, typically center on the use of adaptive equalization and/ 
or dual-branch space diversity.’ 

Most conventional space diversity schemes use either selection 
switching or so-called “in-phase” combining of the diversity branches. 
The latter approach concentrates on maximizing the combiner output 
power rather than on minimizing channel dispersion. Recent work, 


* Authors are employees of AT&T Bell Laboratories. 
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however, has dealt with “out-of-phase” combining, which reduces 
output dispersion by suitably adjusting the relative phases between 
the two branches.*® This approach can completely eliminate dispersion 
for certain two-path propagation situations, but not under more gen- 
eral and realistic conditions. 

In this paper, we describe a combiner in which the relative phases 
and amplitudes of the two branches are controlled. With little increase 
in complexity, this approach allows the effects of both dispersion and 
noise to be jointly minimized. We shall consider this type of combining 
within the context of M-level Quadrature Amplitude Modulation (M- 
QAM) systems. 

Two specific approaches for finding the best amplitude and phase 
adjustments are described in Section II. One in particular, based on 
noncoherent spectrum measurements at the combiner output, is iden- 
tified for further study. The simulation and analysis of this scheme 
are discussed in Section III, and its performance for several postulated 
dual-channel response pairs is assessed in Section IV. Comparisons 
with other diversity and nondiversity approaches are also given and 
attest to the effectiveness of the new scheme. 


ll. THE DIVERSITY COMBINER 
2.1 Rationale 

The idea of using both amplitude and phase adjustments in space 
diversity combining, while not entirely new, has yet to be fully under- 
stood and optimally exploited. We now illustrate the potential power 
of this form of combining under quite general circumstances. To do 
so, we will invoke some recent work on the modeling and analysis of 
multipath fading responses. 

Let H,(f) and H2(f) be the complex frequency responses of a fading 
channel as viewed by two vertically spaced receiver antennas. Under 
nonfading conditions, these functions are flat with frequency at unity 
amplitude (| H,(f)| = | Ho(f)| = 1). In all that follows, we measure f 
from the center of the radio channel, so that H;(f) and H.(f) are 
baseband functions. Moreover, we are interested in their variations 
over just the interval [—W/2, W/2], where W is the channel bandwidth 
in hertz. Since multipath fading arises from a finite number of discrete 
propagation paths, we can present H,(f) and H,(f) in the following 
general forms: 


Ky, 
H,(f) = 2 Ri,exp{—] (wri — O1n)} (1) 
and 


Kg 
A,(f) = Py Ro,exp|—j (wrer — A2x)}. (2) 
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In (1), K; is the number of paths and 7,,, Riz, and 6,, are the time 
delay, amplitude, and phase, respectively, associated with the kth 
path. Similar definitions apply to Ko, Tox, Rox, and 99, in H2(f). 

Typically, multipath propagation on microwave radio links exhibits 
delay spreads on the order of 10 ns or less, i.e., the largest and smallest 
7 values differ by amounts small compared to 1/W, where W < 40 
MHz in the common carrier bands. This observation inspired earlier 
efforts to approximate fading channel responses using low-order pol- 
ynomials in jw.'” 

To be concrete, let H,(/) and H.(f) be represented by the infinite 
power series 


H,(f) = e@*5[A, + jwB, + (Jw)?C, + +++] (3) 


and 
H,(f) = e7**[Ao + jwBe + (jw)?C2 + ---], (4) 


where ¢, and ¢, are arbitrary, and the A’s, B’s, and so on are complex 
coefficients. Using the power series expansion for e’*, they can be 
easily related to the parameters of the functions (1) and (2), e.g., 


Ky . Ky * 
Ay = Dd Rye! B, = ¥ Rye" (ti — 71x), (5) 


and so on. 

The essence of first-order polynomial fitting is this: Given H,(/), a 
value for t; can usually be found such that (A; +jwB;) is the dominant 
part of (3) over [—W/2, W/2], all higher-order terms in Jw being small, 
in some sense. Similarly, a value for te can usually be found that does 
the same for H,.(f), (4). That such first-order polynomial fitting is 
reasonable to do in common carrier channels has been supported by 
theory,’? noncoherently measured data,’ and (more recently) coher- 
ently measured data.’* 

Now suppose that a space diversity combiner were used having an 
adjustable time delay (7) and complex gain (@) in the second branch. 
The composite channel response, as viewed at the combiner output, 
would then be 


H(f) = Hy(f) + Be?*"H2(f). (6) 
If 7 and 8 were adaptively adjusted to be 
7=t,-—t and ~ = —B,/B,, (7) 


we could then write [see (3) and (4)] 


second- and 
H(f) = [A; — A2Bi/Be] + *} higher-order ¢ . (8) 
terms 1n Jw 
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Thus, by proper choice of delay and gain in one branch, the channel 
response could be made dispersionless except for small higher-order 
terms in jw. This is a quite general result for the channels of interest 
and shows the power of the combining approach under study. If 8 
were a phase-only factor [6 = exp(j¢)], such a strong reduction in 
dispersion would only be possible in those fortuitous situations where 
|B. | = | Bel. 


2.2 The combiner structure 


The above discussion suggests both a particular structure, (6), for 
the combiner and particular solutions, (7), for the variable parameters. 
The discussion was intended, however, to provide insight rather than 
to identify a serious design approach. For one thing, a variable time 
delay would be difficult to implement and would offer little incremental 
benefit in most cases. In addition, the solutions of (7) do not properly 
address either the residual dispersion (i.e., higher orders in jw) or the 
receiver noise. If, for example, Ai/B,; and A2/Bz happened by chance 
to be close in value, the first term in (8) would be severely weakened 
in the process of eliminating the jw term; this would enhance the 
effects of both the remaining dispersion and the noise. 

Accordingly, we propose a combiner in which the relative amplitudes 
and phases in the two diversity branches are adaptively adjusted, but 
not the delays. Thus, 


A(f) = BiH (f) + BoH2(f), (9) 


as indicated in Fig. 1, where 8; and @, are adapted gains. (Because the 
dominant thermal noise is introduced before the combiner, several 
ways of adapting 6; and 82 would, in theory, yield equivalent perform- 
ance. For example, one gain could be held fixed, or adapted in ampli- 
tude only, with the other being adapted in both amplitude and phase; 
or both gains could be adapted in both amplitude and phase. Each of 
these approaches would permit the adaptation of the relative complex 
branch gains, which is all that matters.) Moreover, we propose the use 
of control strategies that take proper account of both dispersion and 
noise. 


2.3 Control strategies 


We shall discuss two distinct approaches for controlling the gain 
pair (8,, B2) in (9). The first approximates the theoretically best way 
to do combining when there is no post-combiner equalization. The 
second approach, thought suboptimal, has features that make it at- 
tractive both with and without post-combiner equalization. The new 
scheme reported here incorporates the second approach, whose per- 
formance we will compare with that of the first. 
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Fig. 1—Block diagram of the space diversity combiner. In this scheme, noncoherent 
spectral measurements are made on the combiner output at three or more inband 
frequencies, and the control measure formed from them is used to adjust the complex 
gains 6; and 62. The combining is shown at Radiofrequency (RF) but could be at 
Intermediate Frequency (IF) instead. 


Approach 1: Assume for now that there is no post-combiner equali- 
zation. In that case, the detection bit error rate is virtually minimized 
by choosing (8, 62) to maximize the ratio of sampled signal to root 
mean squared (rms) distortion at the receiver output. By “sampled 
signal,” we mean the half-distance between signal levels, as sampled 
in every period at the in-phase and quadrature detectors; by “rms 
distortion,” we mean the rms sum of thermal noise and Intersymbol 
Interference (ISI) sampled at the detectors. This signal-to-distortion 
ratio [defined formally by (26) and expressed by (28), below] depends 
on H;(f), He(f), 61, 82, and other factors; most important, it is convex 
in 6, and Gy. The result is that there is a unique value of relative 
complex gain, £;/62, that maximizes this ratio. We shall regard as 
optimal any gain pair ((, 2) that exhibits this maximizing relative 
gain. 

There is a practical way to realize optimal values for 6; and Bo, 
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namely, gradient search methods using data decisions.’*’® The action 
is similar to that of an adaptive transversal equalizer, with 6, and B. 
taking the place of the optimized tap gains. We shall refer to this or 
any decision-directed scheme that optimizes ((;, 62), in the sense 
defined above, as Approach 1. 

Approach 2: A potential liability of Approach 1 is that it relies on 
having accurate data decisions, a condition that may not always exist 
(e.g., during recovery from severe fades). The scheme to be reported 
here is based on a different strategy, which we designate as Approach 
2. It consists of (1) performing certain noncoherent spectral measure- 
ments on the combiner output; and (2) sequentially searching over 
(81, B2) so as to maximize a certain quantity [eq. (18), below] computed 
from these measurements. By adapting 6, and 8, in this way, a close 
approximation to the “optimal” condition defined above can be 
achieved, as we show later. That is, Approach 2 should yield near- 
optimal detection performance for receivers having no post-combiner 
equalization. For cases where such equalization is used, Approach 2 
would serve a different purpose—reducing the signal dispersion as 
seen by the equalizer input, thereby simplifying the requirements on 
equalizer design (e.g., number of taps) and improving convergence 
speed. Thus, Approach 2 has the twin virtues of not relying on data 
decisions and having utility both with and without post-combiner 
equalization. 

To describe Approach 2, we refer to the block diagram of Fig. 1. The 
combiner output signal, y(t), has a power spectrum density given by 


Sf) = SCA IAAI? + No(l Bil? + 1 621)| AeA I’, (10) 
Ny 


Signal Noise 


where N, is the power spectrum density of the receiver input noise,* 
Hpr(f) represents whatever receiver selectivity precedes the combiner 
output, and S(f) is the spectral density of the signal (excluding channel 
and combiner effects). More specifically, 


S(f) = So] Hr(f) |7| Hr (f) |’, (11) 


where | H7(f) |? represents the spectral shaping in the transmitter and 
S, is a spectral density scale factor. 

All functions and parameters in (10) and (11) are design-specified 
except | 6,| and | 82], which are controlled by the combiner circuitry, 
and |H(f)|, which must be measured in real time. Our scheme 
estimates | H(f)| at N evenly spaced frequencies (N odd) within the 


* Though not made explicit by the figure, we assume there is sufficient front-end 
amplification that the combiner gains 6, and £2 have no effect on receiver noise figure. 
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channel bandwidth by estimating the corresponding values of S,(/) 
[see eq. (10)]. Based on these estimates, 6, and £2 are adjusted to 
maximize a computed performance measure, Y, which we introduce 
shortly. Before doing so, we define the following: 


Af & Spacing between estimates of | H(f)|, (12) 
where (N — 1)Af s W; 


H 4 Ave{H,,} = “Average Signal Gain”; (14) 
X.ig & Ave{S(f,)}(H)*W = “Signal Power’; (15) 
Xas & Ave{S(f,)(H, — H)*}W = “Distortion Power”; (16) 


Xnoise & No(|8i[* + | B2|")Ave{| Hr(fa) 7} W = “Noise Power”; (17) 


where W and T are the channel bandwidth and digital symbol period, 
respectively. All quantities in these equations are known a priori except 
the H,, values, which are measured. 

We now define the performance measure to be computed and 
maximized, namely, 


Y A pe O.er + X ssise)’ (18) 


This ratio is an approximation, computed from noncoherent spectral 
measurements, of the detector output signal-to-distortion ratio defined 
by (26), below. It is an apt measure to maximize in its own right, for 
the following reasons: In typical digital radio links, noise will not be a 
serious factor unless H,(f) and H.(f) are strongly faded. Therefore, 
maximizing Y will, in most cases, amount to minimizing the ratio 
Xais/Xsig, Which is a measure of the dispersion in H(f). Including Xpoise, 
however, safeguards against minimizing this ratio at an undue cost in 
signal (.X,;,) and thus seriously degrading the signal-to-noise ratio. 
The control strategy is therefore as follows (see Fig. 1): At the 
combiner output, a parallel bank of envelope detectors is used to 
estimate | H(f)]| at N frequencies. The spectral samples are digitized 
and applied to a microprocessor, which computes Y. This measure 
drives the search over |{;| (or |@2|) and @ = Arg{@2}, i.e., these 
quantities are adjusted so as to maximize Y. Typically, they are 
adjusted iteratively, e.g., ¢ is changed in 0.1-radian step until a local 
maximum is found; then | 6,| or | 82| is changed from 1 in steps of 0.1 
until a maximum is found; and this process repeats, possibly using 
smaller steps in successive rounds, until Y can no longer be increased 
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by varying either 6 or ¢. If each measurement (i.e., set of estimates of 
| H(f)| at N frequencies) and computation for Y takes ty seconds, 
and N, steps are needed to find (f;, 82), then the “solution time” of 
the combiner will be about N,ty. This number should be small 
compared to 1 second to achieve timely adaptation to multipath fades. 


2.4 Measurements 


We now discuss the scheme for measuring the set of H,,’s in (13). In 
our simulations (Section III) we treat only the case where N = 3 and 
Af = W/2 (i.e., three samples, taken at the channel edges and center). 
In practice, the outer samples would probably be closer-in so as to 
minimize errors from adjacent-channel interference. Also, higher val- 
ues of N (e.g., N = 5) might be worthwhile. 

To see how accurate estimates of H, might be obtained, let G(/) 
represent a low-pass power gain function with bandwidth b/2 «< W. 
We envision the measurement of H, as involving a bandpass filter 
with power response G(f — nAf) followed by envelope detection and 
ty-second averaging of the detector output (Fig. 1). Referring to eq. 
(10), the average power at the output of the bandpass filter will be 


F, = f IS(PIH(AL 
+ NABI? + |B2I Ha IIG(f— nAfdf, (19) 


where S(f) is defined by (11). Assuming a square-law detector, the 
time-averaged detector output will be 


(P,) =P, + | Fluctuation ea 


: | Variance ~P,,/bty (20) 


We now define two constants related to the system design functions, 
namely, 


tn & No i) | Hr(f) ?G(f — nAf)df (21) 
and 


i, | S(AG(f — nAfodf. (22) 


Based on (19) and the fact that H(/f) changes little over the bandwidth 
b, a microprocessor can estimate H,, using the formula 


a, = \ / Pad = toll Bl? + 18217) 


n 23 
e (23) 
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where (P,,) is measured in real time; | 61] and | $2| are controlled 
parameters of known value; and 7, and ¢, are predetermined constants. 

Equation (23) shows how A, is computed in terms of measured or 
known quantities. To see what this computed number represents, we 
insert (19) through (22) into (23) and obtain 





{ SING RAMIHPPof 7 
j= [— {Rime stey ly 
{ S(PG(f — nApyaf fluctuation noise 


We can now cite choices for b and ty that lead to accurate and 
sufficiently rapid estimations of H,,. As b gets very small, the first 
term under the radical sign in (24) approaches H2, so that the major 
inaccuracy in H,, is due to fluctuation noise. To be more precise, the 
first term is close to H2 so long as | H(f) |? changes little over the 
passband of G(f — nAf). Since we are considering propagation media 
with delay spreads of just a few nanoseconds, the design rule b < 2 
MHz should permit more than adequate resolution in this regard. To 
achieve low mean-square fluctuation noise as well [second term in 
(24)], the condition bty = 4000 should be satisfied [see (20)]. Thus, 
with b = 2 MHz and ty = 2 ms, H,, can be approximated with high 
accuracy by the quantity H,. Moreover, this design choice would 
permit numerous iterations of the search over £; and #2 before the 
medium response changes appreciably. 


lil. PERFORMANCE STUDY 
3.1 General 

We have written a set of computer programs to simulate the behavior 
of the combiner scheme described above (Approach 2) and to analyze 
its performance and that of other receiver techniques. Each simulation 
is done for a specific pair of fading functions, H,(f) and H2(/f), and for 
a specific value of Carrier-to-Noise Ratio (CNR). What is simulated 
is the sequential search over 6; and $2, as performed by a receiver in 
real time to maximize the computed measure Y [see eq. (18)]. 

The analysis programs compute a detection performance measure 
for a nonequalized receiver using the (8, 6.) pairs derived in the 
‘simulations for Approach 2. The same measure is also computed for 
the cases of optimal combining (Approach 1) and no combining (non- 
diversity). Also, the analysis programs examine the signal dispersion 
at the combiner output for these various cases, which has relevance to 
receivers with post-combiner equalization. 


3.2 Response pairs studied 


We have specified eight distinct pairs of H,(f/) and H2(f) for pur- 
poses of study. These pairs are collectively representative of what 
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Table |—Dual-channel response pairs studied 


PLOTS OF PATH Ty, 
lH, (A), 142(f)| |DIVERSITY| NUMBER, IN ns) R,expO,) 
(IN dB) k 
4 HIN dB 
-4.0 1.0 + j0.1 
+1.0 -0.5 + j0.9 
-2.0 1.1 +/0.4 
+4.0 0.9 - j0.4 
0.2 + j0.1 
~5.0 -1.1 + (0.6 
+3.0 0.5 + /1.1 


1.0 + /0.1 
-0.9 + j0.5 


1.1 + 0.9 


-1.0 + j0.5 


1.0 + j0.0 
-0.9 + (0.0 


1.1 + 0.0 


-0.5 + j0.0 





might arise in actual radio links using space diversity. In each case, 
H,(f) corresponds to a one-, two- or three-path medium, and similarly 
for H.(f). [In terms of (1) and (2), K,; = 1, 2, or 3 in each case, and 
similarly for K2.| The corresponding time delays and complex gains 
are summarized in Table I. Also shown for each case are graphs of 
| H,(f)| and | H2(f) |, in decibels, over a 40-MHz bandwidth. 
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Table I—(Cont.) Dual-channel response pairs studied 


PLOTS OF SPACE PATH T; 
lA (Fl, lHa(f)l [| DIVERSITY| NUMBER, (IN ns) R,exp(/§,) 
(IN dB) CHANNEL k 
-2.0 1.0 + {0.0 
+2.0 -0.872 - (0.223 
-1.0 0.54 + /0.84 
+3.0 -0.66 - j0.60 


1.0 + (0.0 


-0.5 + j0.0 


~3.0 -0.068 - (0.998 
+1.0 0.48 + /0.76 


0.0 - /1.0 
0.5 + j0.5 
1.0 + (0.0 


0.0 - /1.0 
-0.5 + j0.5 
0.5 + /0.0 


-1,0 + j0.5 
0.2 + 0.7 
0.5 - j1.0 


1.0 - j0.2 
1.0 - /0.2 
0.5 + /0.5 





Because of additive noise, the receiver performance for a given 
response pair would be affected by any amplitude scaling of (H,(/), 
H.(f)). We permit the possibility of such a scaling, for each case in 
Table I, by including it in the carrier-to-noise ratio parameter dis- 
cussed below. 


SPACE DIVERSITY COMBINING 895 


3.3 Signal design assumed 


We assume throughout the study that the signal spectral density 
S(f), (11), is rectangular over the channel bandwidth, with magnitude 
S,. This corresponds to the use of ideal Nyquist pulses with 0-per- 
cent roll-off factor (a = 0) and a symbol rate equal to the bandwidth 
(1/T = W). While not realistic, this assumption both simplifies the 
analysis and leads to somewhat poorer results than for other (more 
tapered) spectra.’” It thus serves both convenience and our confidence 
in whatever favorable outcome the study predicts. Adding to the 
“worst-case” nature of the results is that they are obtained for the 
largest common carrier channel bandwidth, W = 40 MHz. 


3.4 Carrier-to-noise ratio 


Finally, we introduce a so-called flat fading gain, g, with which we 
can amplitude-scale the various pairs (H;(f), Ho(f)) in Table I. Ac- 
cordingly, we define the flat-fading carrier-to-noise ratio to be 


CNR 4 g*S,/No. (25) 


A typical system value for this quantity, with g? = 1, is 10° (60 dB). 
We will consider a wide range of values for CNR, thereby accommo- 
dating a wide range of amplitude scaling factors. 


3.5 Method of analysis 


We consider an M-level QAM system with the previously noted 
Nyquist signaling. At each baseband detector, data decisions are made 
every T' seconds by comparing the sampled input with a set of decision 
thresholds. We now introduce the following new quantities: 

Ps & The squared signal sample (excluding ISI and noise) at 
either baseband detector, when a data value of +1 or —1 is 
being detected; 

P;, Py 2The mean-squared ISI and noise, respectively, associated 
with the periodic samples at either baseband detector. 


We define the detection signal-to-distortion ratio to be 
pp & Ps/(P; + Py) (26) 


and note that this quantity yields a reasonably tight upperbound to 
the bit error rate, via’® 


Bit Error Rate < 2 exp(—pp/2). (27) 


We thus regard pp as a proper index for evaluating digital radio receiver 
performance. 

Now consider a space diversity receiver using no post-combiner 
equalization and exhibiting a channel response H(f), (9), at the 
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combiner output. For the assumed M-QAM system, with optimal 
timing and carrier recovery, we can show that 

















3 CNR(H72.,x) ee: r, (28) 
Pop = ee ee Oe, ee , 
M—1CNR(|A[? — Hox) + (161? + [62/2) “@-1 
where 
Ww/2 
Hinax & Max| T H(f)e?°df (29) 
t, -w/2 
and 
W/2 
[HF 4 r{ | H(A) Pf. (30) 
~Ww/2 


Note that the dependence of pp on M resides entirely in the factor 
3/(M — 1). By using I as our signal-to-distortion measure, therefore, 
we can remove M from the set of problem variables. 

Since H(f) is linear in 8; and $, a solution exists for these gains 
(actually, for their ratio) that maximizes I’. The solution can be found 
analytically by a variation on the method outlined in Section 5.1 of 
Ref. 15. In an actual receiver, the solution can be closely realized using 
practical circuitry, which we have designated as Approach 1 (Section 
2.3). This solution for (6;, G2) is defined here to be optimal, and we 
will compute and present I results corresponding to it. These will be 
compared with I results for the “suboptimal” (61, 62), which we define 
as the pair produced under Approach 2 (Section 2.3) and which are 
obtained here via computer simulations. For completeness, we will 
also show I results for the case of no diversity, i.e., either 6; = 1, 
B82 = 0 (Branch 1 only is processed) or 6; = 0, 8. = 1 (Branch 2 only is 
processed). 

Finally, in consideration of receivers using post-combiner equaliza- 
tion, we define a “dispersion index” for H(f) as follows: 


I> 


fo H(f) |} 
R f 


SSS <= 

sera | < 20 MHz. (31) 
This is just the range of | H(f)| over a 40-MHz channel bandwidth. 
Gersho has shown that adaptive equalizers converge more rapidly 
when this ratio is close to unity.’? Another likely benefit of near-unity 
R is that effective equalization should be attainable using a relatively 
small number of equalizer taps. We will compute R for the same cases 
that we compute I. 
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3.6 Simulation of the search algorithm for Approach 2 


The simulations assume three spectrum measurements located at 
—20 MHz, 0 MHz, and +20 MHz relative to the band center. As noted 
before (Section 2.4), we can identify design parameters that drive 
measurement inaccuracies to negligible levels. For this reason, we 
assume exact spectral estimates and have not attempted to simulate 
measurement inaccuracies. 

For a given response pair (‘Table I) and CNR, (24), the simulation 
proceeds as follows: During the search process one of the 8 values is 
always set to unity with the other value less than or equal to unity. 
Initially, | 8:,||@2| and @ are set to 1.0, 1.0 and 0, respectively. Then 
@ is varied in 6¢-radian steps until Y, (18), is maximized. Next, | G2 | 
is varied in steps of 66 (if the simulation calls for increasing | 82 | above 
unity, we set |@.| = 1 and decrease || instead) until Y is again 
maximized. The procedure is then repeated until Y cannot be further 
maximized. 

In the simulations, we generally started with 68 = 0.1 and 6¢ = 0.1 
radian. We noticed that a steady-state result usually required about 
twice the number of steps needed to reach it directly. For example, if 
the solution were | 8; | = 1.0, | B2| = 0.4, and ¢ = 3.1 radians, a “direct” 
path would entail six 0.1-step changes in |8.| and 31 0.1-radian 
changes in @, or 37 steps in all. The maximizing algorithm described 
above, however, was found to take about 74 steps. Using 68 = 0.01 and 
5@ = 0.01 radian would, of course, require 740 steps. Our approach 
was to use 0.1 for both increments until a stable solution was reached 
and then to proceed to a finer solution by changing both increments 
to 0.01. The end results using this faster two-stage process were 
identical to those obtained by using 0.01 throughout. 

Using the two-stage process, the measurement steps needed to reach 
a solution should not exceed perhaps 200 for a “cold start” adjustment 
of the combiner gains. For ty = 2 ms per measurement, this corre- 
sponds to a solution time of 0.4 second. In the dynamic situation 
where the receiver tracks the channel variations, the solution time 
would be much shorter. 

An important question is whether the derived solutions for (81, Be) 
are unique. We explored this question for each of the eight cases in 
Table I, plus a few others. For each case, we simulated the search 
process using 24 different starting points for | 8,|, | G2| and ¢. Specif- 
ically, (| 61], | 82|) was initialized to each of three different pairs of 
values, namely, (0.5, 1.0), (1.0, 1.0), and (1.0, 0.5); and, for each such 
pair, ¢ was initialized to each of eight values, namely, 0, 45, 90, --- 
315 degrees. In studying the 24 solutions for each case, we observed: 

1. The solution for (6;, G2) was not unique in every case; for some 
cases, multiple local maxima for Y were found to exist, and the search 
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algorithm converged to one or another, depending on the starting 
point. 

2. The multiple-solution condition appears to exist only at high 
CNR (240 dB); for CNR = 20 GB, all solutions for (61, G2) were found 
to be unique.* 

3. In those cases where multiple solutions did exist, the lowest of 
the resulting Y values was always quite good (>30 dB) and usually not 
much below the highest Y. . 

From these observations, we tentatively conclude that the proposed 
scheme converges to stable solutions under all circumstances. Further, 
the solutions seem to be nonunique only in circumstances where all of 
the possible solutions are satisfactory. 


IV. RESULTS 


A qualitative assessment of the combining scheme can be gained 
using graphical results. Figures 2, 3, and 4 show, for each of Cases 1, 
3, and 4 in Table I, plots of Hi(f), He(f), and the H(f) obtained when 
CNR = ~. In each figure, the severe dispersions of H,(f) and H.(f) 
are seen to be almost eliminated by the combining: H(f) presents a 
nearly flat amplitude response and a nearly linear phase response over 
the channel bandwidth. 

To make the assessments more quantitative, we have computed I 
and R, (28) and (31), for each case in Table I; for each of the CNR 
values 60 dB, 40 dB, and 20 dB; and for each of the following four 
solution strategies for (81, 82): 

Solution 1: Optimal; based on Approach 1, it leads to a maximum 
for I with respect to 6, and Bo. 

Solution 2: Suboptimal; based on Approach 2, it leads to a maximum 
for Y with respect to 6, and fo." 

Solution 3: Nondiversity; Branch 1 processed only. 

Solution 4: Nondiversity; Branch 2 processed only. 

We define I; and R;(i = 1, 4) to be the decibel values of I and R, 
respectively, for Solution 1. The results are given by Tables I, III, 
and IV, each corresponding to one of the three values of CNR. 


V. DISCUSSION AND CONCLUSION 


The data of Tables II through IV lead us to the following observa- 
tions: 
1. Comparing I, with P'3 and Ty reveals that an unequalized receiver 


* Under low CNR conditions, the noise term dominates. Since the noise term is 
quadratic in £,; and fz, it leads to a unique maximum for Y. 

TIn all cases where multiple solutions were found to exist, we used the one corre- 
sponding to the higher Y. 
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Fig. 2—Relative amplitude and phase plots for Hi(f) and H,(f), corresponding to Case 1 in Table I, and the resulting combiner output 
response, H(f), when (; and #2 are derived using the control scheme in Fig. 1; no noise present. The plots cover a 40-MHz bandwidth and, 
for convenience, all magnitudes and phases are shifted so as to be zero at band center. 
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Fig. 3—Same as Fig. 2, for Case 3 in Table I. 
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Fig. 4—Same as Fig. 2, for Case 4 in Table I. 


Table Il—Combiner performance results for CNR = 60 dB 


Detection Measure (in dB) Decibel Range of |H(f)| 
Case Tr; Tr, I's r 4 R, Re Rs Ry 
1 43.80 42.67 11.12 20.88 0.11 0.062 9.34 2.77 
2 46.99 46.99 46.99 11.15 0.006 0.0 0.0 9.79 
3 36.06 35.62 9.17 13.39 0.21 0.14 12.25 6.82 
4 53.13 53.12 6.82 25.91 0.002 0.001 11.35 1.39 
5 29.23 27.60 2.80 3.02 0.37 0.27 16.96 16.49 
6 35.32 34.99 18.21 5.38 0.26 0.21 3.19 19.09 
7 42.28 40.72 18.31 17.80 0.11 0.13 3.76 3.08 
8 63.84 59.24 10.70 29.02 0.006 0.022 7.38 1.06 


Table I!]—Combiner performance results for CNR = 40 dB 


Detection Measure (in dB) Decibel Range of |H(f)| 
Case rT 1 IT. T's T, R, R, R3 R, 
1 .41.01 40.48 11.12 20.86 0.11 0.06 9.34 2.77 
2 27.46 26.99 26.99 11.15 0.44 0.0 0.0 9.79 
3 29.91 29.60 9.16 13.39 0.41 0.18 12.25 6.82 
4 33.388 33.33 6.77 25.51 0.11 0.05 11.35 1.39 
5 27.11 25.70 2.79 3.00 0.40 0.27 16.96 16.49 
6 31.29 30.97 18.13 5.37 0.36 0.21 3.19 19.09 
7 38.08 37.72 18.30 17.71 0.17 0.13 3.76 3.08 
8 46.23 46.00 1068 28.96 0.022 0.014 7.38 1.06 


Table IV—Combiner performance results for CNR = 20 dB 


Detection Measure (in dB) Decibel Range of |H(/)| 
Case Yr, T, T3 YT, R, R, Rs R, 
1 25.98 25.76 10.68 19.67 0.58 0.43 9.34 21 
2 12.91 12.54 6.99 10.96 5.72 4,31 0.0 9.79 
3 14.68 14.35 8.29 12.93 3.96 2.85 12.25 6.82 
4 15.67 15.63 3.58 15.64 1.54 1.37 11.35 1.39 
5 11.72 11.68 1.97 2.19 2.16 1.83 16.96 16.49 
6 15.21 15.14 13.49 4.76 1.99 1.62 3.19 19.09 
i 23.80 23.74 17.69 12.75 0.75 0.63 3.76 3.08 
8 27.63 27.63 9.09 25.41 0.22 0.20 7.38 1.06 


with “suboptimal” space diversity combining would generally perform 
much better than one using no diversity (13 or 4) or selection diversity 
(the larger of I's and T'4). 

2. Comparing I, with TF, reveals that the “suboptimal” combiner 
would perform nearly as well as the “optimal” one, at least for the 
response pairs considered here. 

3. For the “suboptimal” combiner, the range (R2) of the post- 
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combiner frequency response is generally much smaller than that for 
either branch alone (R3 and R,). It is also comparable to that for 
“optimal” combining (R,), often being even smaller. 

4. Comparing Tables II, III, and IV reveals the trends of combiner 
performance as CNR decreases from high values (=40 dB) to low ones 
(20 dB). At high values, the combiner emphasizes minimum dispersion, 
as reflected in the data for R; at low values, combiner action aims 
more at minimizing noise and so dispersion reduction is limited. 

It thus appears that the noncoherent measurement/control scheme 
described here would be effective in any receiver situation. This 
includes receivers without post-combiner equalization, wherein I 
should be maximal to optimize detection; and receivers with equali- 
zation, wherein R should be minimal to facilitate equalizer conver- 
gence. It is clear that both aims are served by the same algorithm. 
Moreover, the scheme is simple, fast-acting, and operates in proximity 
to the combining circuitry. 

Some issues remain to be settled. One is whether the search algo- 
rithm will always converge to either the near-optimal solution or to 
one that, in any case, yields very good performance. The results 
obtained for our limited sampling of response pairs offer encourage- 
ment on this score. Another issue concerns the effectiveness of diver- 
sity combining over a statistical ensemble of fading conditions, that is, 
we have obtained results here for just eight selected response pairs. In 
a separate study, however, we invoked a recently developed statistical 
model for dual-diversity channels”°** to simulate a large population 
(ensemble) of response pairs (H,(f), Ho(f)), and we obtained proba- 
bility distributions for I and R over that population. The results are 
reported in a companion paper.” 
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A Simulation Study of Space Diversity and 
Adaptive Equalization in Microwave Digital 
Radio 
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In this paper we analyze the performance of M-level quadrature amplitude 
modulation digital radio systems subjected to microwave multipath fading. 
We consider two kinds of adaptive receiver techniques, either singly or in 
combination: dual space diversity and adaptive equalization. The space diver- 
sity is assumed to be of either the selection type or the continuous-combining 
type, and the equalization is assumed to be ideal. We describe a specific form 
of combining which is optimal when no post-combiner equalization is used. A 
primary aim of the study is to quantify the performance of this combining 
approach and to compare it with alternate strategies. The study uses Monte 
Carlo simulations of the dual-channel fading response functions based on a 
recently published statistical model. For each response pair generated, a 
receiver detection measure is derived analytically in terms of the system 
parameters and receiver approach. Probability distributions of this measure, 
obtained by simulating several thousand response pairs, are then computed. 
They can be interpreted as displaying the link outage probability as a function 
of the number of modulation levels (M). We find that the appropriate 
combining scheme can serve in some cases to avoid the need for adaptive 
equalization. Also, where post-combiner equalization is used, the same scheme, 
while no longer optimal, can sharply reduce the dispersion seen by the equalizer 
input. 


I. INTRODUCTION 


A continuing challenge in microwave digital radio is to find ways to 
counter multipath fading. Many recent efforts—theoretical, experi- 


* Authors are employees of AT&T Bell Laboratories. 
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mental, and developmental—have addressed the use of space diversity 
and adaptive equalization, individually or in combination, to main- 
tain stringent outage objectives for increasingly high-level modula- 
tions.’"* This paper reports on a comparative study of several receiver 
processing approaches, enlarging on the evaluations presented in a 
companion paper.” 

The receiver processing approaches studied here consist of three 
space diversity options, namely, no diversity, selection diversity, and 
continuous-combining (or combining) diversity. For each of these op- 
tions, we evaluate performance statistics both with and without the 
presence of ideal adaptive equalization. 

Our investigation is facilitated by a new statistical model for the 
response functions of two space diversity channels.’°’” We describe 
that model in Section II and introduce some extensions germane to 
the present study. 

We assume throughout this work a system using M-level Quadrature 
Amplitude Modulation (M-QAM) with cosine roll-off spectral shaping. 
This system and the various space diversity/equalization configura- 
tions to be studied are described in Section III. 

For any given dual-channel response pair and receiver processing 
approach, we can compute one or more relevant performance meas- 
ures. By combining the published model with Monte Carlo simulation 
techniques to generate response pairs, we can obtain probability dis- 
tributions of these measures over an ensemble of multipath fades. In 
Section IV we define the performance measures to be studied and 
outline the analytical methods by which they are computed, and in 
Section V we describe the method of Monte Carlo simulation and 
present quantitative results. These are in the form of probability 
distributions and reveal the influence of various system parameters, 
channel parameters, and receiver processing approaches on link per- 
formance. 


Il. DUAL-CHANNEL FADING MODEL 
2.1 Background 

We are dealing with propagation media and systems wherein the 
multipath delay spread generally is small compared to the inverse 
channel bandwidth. Consequently, fading responses vary smoothly 
with frequency over a given channel and can be accurately approxi- 
mated using simple functions. A popular complex function for this 
purpose is the one introduced by Rummler, initially for nondiversity 
channels’®*!? and, more recently, for dual-diversity channels.’*” In its 
most general form, Rummler’s function approximates a complex re- 
sponse H(f) on a given frequency interval, [-W/2, W/2], with W < 
40 MHz, by 
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H(f) = [e, + coe I" Je eto, lfl <= W/2, (1)* 


where c,; and cz are complex coefficients that vary slowly with time 
but can be regarded as quasistatic; 7 is a fixed parameter (6.3 ns); and 
t, is a time delay (nominally, the propagation delay through the 
medium for which H(f) is the response). If H(/) were precisely known 
on [—W/2, W/2], H(f) could be fitted to it by first choosing least- 
mean-squared values for c; and Cc, (i.e., by minimizing, for given ¢,, the 
integration of | H(f) — H(f)|? over [—W/2, W/2]), and then minimiz- 
ing the result with respect to f,. 

An alternate form for H(f) that more closely resembles the form 
used by Rummler is 


H(f) = a[1 — beW" Joie, | f| <= W/2, (2) 


where a is real and positive; b is complex; | b| < 1; and 7 is either +6.3 
ns (“minimum-phase” response) or —6.3 ns (“nonminimum-phase” 
response). For each of these two conditions, it is a simple matter to 
relate a, b, ®,, and t/ to ci, co, and ¢, in (1). The question “minimum 
phase or nonminimum phase” has been left open to date, since only 
amplitude versus frequency data have informed the major attempts at 
statistical modeling.** ”° . 

In Rummler’s work, the phase factor exp(j(®, — wt3)) is omitted. 
Again, this is necessitated by the absence of phase versus frequency 
information in the available databases. In nondiversity reception, 
moreover, this factor is immaterial since the phase ©, and time delay 
t, would be tracked, and their effects removed, by the carrier and 
timing recovery circuits of the receiver. 

Now, however, consider the dual-diversity case. As before, each of 
the two responses, H,(f) and H.(f), can be approximated on [—W/2, 
W/2] using functions like (1) or (2). Since combining diversity involves 
forming the composite response §,Hi(f) + §8.H2(f), the ®, term for 
each channel is immaterial; these phases can be regarded as absorbed 
into the phases of 8; and Bo. The difference in t for the two channels, 
however, is another matter. We will return to this point later, after 
we summarize the statistical model. 


2.2 Joint statistics: Rummler’s model 


The dual-diversity model is derived from data collected on a 26.4- 
mile path in Georgia in the 6-GHz band. The details—data, methods, 
and results—are well documented in Refs. 16 and 17. Here we will 
merely state some major results, using slightly different symbols where 
appropriate. 


* Throughout this work, f is frequency measured from the center of the RF or IF 
channel under study, so that all responses are equivalent low-pass functions. 
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To begin with, the expected number of fading seconds on a micro- 
wave hop per year is 


T, = 52,800c(F/6)(D/25)°, (3) 


where F is the microwave carrier frequency in GHz; D is the hop 
length in miles; and c is the terrain factor, varying between 0.25 and 
4, During the rest of the year (normal propagation), the two response 
functions are H,(f) = H2(f) = 1 + j0. During multipath fading, 
however, they are 


Ay(f) = a[1 — byexp(+jw7)], (4a) 
H2(f) = a2[1 — brexp(+jwr)], (4b) 


where 7 = 6.3 ns and the proper sign to use before it is an open 
question. Both functions apply over the limited range |f| < W/2, 
where W < 40 MHz. 

The joint statistics of the a’s and b’s in (4) were published initially 
in Ref. 16. To smooth out some apparent artifacts of the data collec- 
tions and reductions, Rummler subsequently published a “rational- 
ized” version in Ref. 17. We will not reproduce his mathematical 
descriptions for the joint statistics, but a few important features should 
be noted: 

e The four quantities | b,|, Arg{b,}, |b2| and Arg{b.} are mutually 
independent random variables, which means that the relative 
shapes of H,(f) and H.(f) are statistically independent. 

e The amplitude factors a; and az are lognormal random variables, 
1.e., their decibel values are Gaussian. Moreover, the mean of each 
Gaussian variable is a function of the magnitude of the corre- 
sponding 0; each of the variations about the mean has a standard 
deviation near 7.0 dB; and the two variations are correlated. 

e For the data reductions reported in Ref. 16, the correlation factor 
(po) between these variations was 0.65. In Ref. 17, Rummler gen- 
eralizes this result by giving a simple empirical relationship be- 
tween p and the vertical spacing of the diversity antennas. 


2.3 Extensions and sensitivity considerations 


We have modified Rummler’s model, for purposes of our study, in 
three ways. First, we address the minimum-phase/nonminimum-phase 
question in the following manner: Let 


_ Ja[1 — b\e%*"] with probability p 
Hy(f) = a — e/*"] with probability (1 — p). (5) 


Assuming that 7 = 6.3 ns and |b| < 1, we see that the first form is 
minimum phase and the second form is nonminimum phase. Thus, p 


910 TECHNICAL JOURNAL, APRIL 1985 


is the fraction of fades for which the form of H;(f) is minimum phase. 
We use the same description for H2(f), and assume that the minimum- 
phase condition occurs independently for the two response functions. 
In our simulations, we treat p as a parameter and vary it from 0 to 1 
to assess its importance. 

Second, we assume that H.(f) contains an additional phase factor, 
exp(—jw6t), where 6¢ is the difference in t for the “best” approxima- 
tions to the two channel responses. This quantity cannot be known 
without simultaneous coherent measurements on both channels. We 
estimate, however, that it can be +2 ns just from time delay misalign- 
ments between the two diversity branches. Allowing as much as +4 ns 
for propagation differences, we speculate that 6t lies between +6 ns. 
We parameterize it accordingly in our simulations, always holding it 
fixed for a given run of Monte Carlo trials. 

Finally, we allow the possibility that the correlation factor (p) 
between the decibel variations of a, and a, can differ from one path or 
time period to another, even for the same antenna spacing. We treat 
it parametrically in our simulations, varying it from 0 to 1 to assess 
its importance. 

To summarize, we are concerned that our simulation results be 
“robust”, i.e., not sensitive to—or artifacts of—uncertain features of 
the model. For this reason, we have identified some key uncertainties, 
cast them in terms of numerical quantities (p, 6t, and p), and treated 
these quantities as parameters. As for the joint statistics assumed for 
the a’s and b’s in (4),'*"” it is conjectural whether these descriptions 
would be applicable to all paths and time periods. For now, we regard 
these statistics as representative and useful for study purposes. 


Hil. RECEIVER PROCESSING APPROACHES 
3.1 General receiver structure 


All the receiver processing approaches to be considered here can be 
cast as special cases of the general arrangement in Fig. 1. The original 
transmission is an M-QAM signal whose baseband pulse has the 
Fourier transform VC(f), with C(/) being the cosine roll-off func- 
tion.” The two receiver inputs differ only in that the propagation 
medium delivers a different response into each space diversity antenna. 
Typically, the antennas would be vertically displaced by 40 feet or 
more so as to receive relatively uncorrelated responses during fading. 
Each diversity branch is shown having a separate, adaptable complex 
gain. The space diversity combiner is terminated by a signal adder 
that delivers a single-branch signal to the remainder of the receiver. 
Although shown here at RF, this circuitry could as readily be at IF. 
Also (though not made explicit by the figure), we assume there is 
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Fig. 1—Block diagram of a digital radio receiver with dual-channel space diversity. The modulation is M-QAM and the spectral shaping is cosine 
roll-off. 


sufficient front-end amplification that the combiner gains 6; and £2 
have no effect on receiver noise figure. 

The combiner output is shown applied to an adaptive equalizer, with 
response H.,(/). For the cases studied that do not involve the use of 
equalization, we will set H.,(f) to unity. The second post-combiner 
filter has a fixed response, VC(f), which provides matched filter 
reception in the absence of fading. Because the receiver is linear, the 
order of the adaptive and fixed filters is immaterial, as is the manner 
in which each is distributed among RF, IF, and baseband stages. 

The coherent demodulation stage uses a quadrature pair of local 
oscillator signals which we assume to be optimally phased.* The 
baseband outputs are sampled in the detector every T' seconds, where 
T is the symbol period (1/T is the baud), and we assume the sampling 
phase (timing epoch) to be optimal. 

What is not shown in Fig. 1 is the means for controlling the 
adaptable gains, i.e., the tap gains of the equalizer (if used) and the 
diversity branch gains, 8; and 62. Regarding the former, adaptive 
equalization is a well-developed art, and so we shall merely assume 
idealized equalizer responses, confident that they can be closely real- 
ized in practice. Regarding the latter, we shall specify idealized criteria 
for determining £; and (2, both for selection diversity and combining 
diversity, and shall apply these criteria analytically in our performance 
calculations. Only in the case of combining diversity, however, will we 
address the issue of practical implementation. 


3.2 Space diversity options 


We delineate the three particular diversity options as follows: 

1. Nondiversity, wherein only the response H;(f) (corresponding to 
the higher antenna) is received; thus, (81, 82) = (1, 0). 

2. Selection diversity, wherein (8;, 8) is either (1, 0) or (0, 1), 1.e., 
only one branch is processed at a given time, depending on which is 
deemed “best” in some sense. 

3. Combining diversity, wherein 8, and $2 are continuously variable 
and are adjusted according to a specified criterion. 

We will discuss these options first for the case of no adaptive equali- 
zation (H.,(/) = 1 +0 in Fig. 1). Accordingly, the composite channel 
response as seen at the fixed filter input is 


H(f) = Biff) + BeA2(f), (6) 


where all quantities are complex. 
The nondiversity case is included here for purposes of comparison. 


* Optimal carrier and timing phase is meant, in these discussions, to denote phases 
for which the detection signal-to-distortion ratio is maximal. 
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It will be shown to yield disastrous performance results in the absence 
of equalization, as expected. 

For the case of selection diversity, the signal adder would, in practice, 
be a selection switch operated by real-time decisions of some kind. We 
will obtain results for a theoretically optimal selection strategy, de- 
scribed as follows: At all times, (6;, 62) is either (1, 0) or (0, 1), 
whichever yields the lower bit error rate for the prevailing response 
pair, (Hi(f), H2(f)). Since selection diversity is included here primarily 
for purposes of comparison, we ignore the issue of practical implemen- 
tation. 

The combining diversity case is the one that interests us most. We 
will show that combining diversity alone (i.e., without equalization) 
can go a long way towards combatting dispersion as well as noise. As 
explained in Ref. 15, this potential exists for any channel over whose 
bandwidth the frequency response can be approximated by a first- 
order polynomial in jw. This condition applies to the microwave 
common-carrier channels at hand.” 

The criterion we will assume for adapting (6, 82) in the absence of 
equalization is maximization of a particular signal-to-distortion ratio 
at the detector, where “distortion” includes Intersymbol Interference 
(ISI) and noise. This ratio, which we call pp, is defined and analyzed 
in Section 4.1. We will see that pp is quadratic in 8, and £2; thus, a 
unique maximum exists and the (fi, 62) pair that produces it can be 
found analytically. More important, it is well known that this solution 
for (81, B2) can be realized in a receiver by means of practical circuitry, 
specifically, via decision-directed gradient search algorithms.”””° 

In Ref. 15 we describe an alternative scheme for adapting ((;, Be) 
that does not rely on data decisions and yet produces results close to 
those using the more conventional approach. The scheme involves 
measuring the combiner output power spectrum density at three or 
more inband frequencies, computing from these samples a measure 
that approximates pp, and using this measure to drive the search over 
By and Bo. 

With all this in mind, we will derive analytical solutions for the 
(81, B2) pair that maximizes pp in the absence of equalization, and will 
denote it by (81, B2). Since maximizing pp minimizes an upper bound 
on the bit error rate (Section 4.1), we will call this gain pair optimal. 
Finally, we will assume that optimal or near-optimal solutions can 
always be realized in practical receivers using schemes such as those 
mentioned above. 


3.3 Receivers with adaptive equalizers 


A variety of practical approaches could be assumed for the equalizer 
in Fig. 1. These include the reciprocal equalizer,”4 the Minimum Mean 
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Square Error (MMSE) equalizer,”’ the Decision Feedback Equalizer 
(DFE),”° and adaptive cancellation.”’ In terms of detection signal-to- 
distortion ratio, these four approaches improve in the order cited. 
However, in terms of outage performance over a large ensemble of 
fading responses, there seems to be little to differentiate among 
them.” We shall therefore assume an ideal reciprocal equalizer, both 
because it is the easiest to analyze and because it provides a tight 
worst-case bound on the performance of the more optimal equalizers. 
The reciprocal equalizer has a response © 


1 1 
H(f) BiHa(f) + BoHo(f)’ 


Its obvious effect is to restore the received signal to what it would 
have been without multipath fading and space diversity. In the process, 
it eliminates ISI but produces a noise enhancement proportional to 
the integral of C(f)/| H(f) |° (see Fig. 1), where C(f) is the cosine roll- 
off function. We will analyze this equalizer in conjunction with the 
three space diversity options (nondiversity, selection diversity, and 
combining diversity) discussed in Section 3.2. 

In evaluating selection diversity with equalization, we will again 
assume a theoretically optimum strategy: At all times, (81, 82) = (1, 0) 
or (0, 1), whichever yields the lower bit error rate for the prevailing 
response pair, (H;(f), H2(f)). As before, we will not consider the issue 
of practical implementation. 

In evaluating combining diversity with equalization, we will consider 
three practical strategies for adapting 8, and G2. They are as follows: 

1. The relative amplitude and phase of 8, and £8, are adjusted to 
maximize the ratio of signal power to noise power at the combiner 
output. 

2. Only the relative phase of 8, and 2 is adjusted to maximize the 
same ratio, with | @;| = | 6.| = 1. 

3. 8, and f» are adjusted to those values, 8, and {, that would 
maximize pp in the absence of equalization. 
Missing from this list of strategies is the (G;, G2) pair that maximizes | 
pp in the presence of equalization. Unfortunately, the analytical solu- 
tion for this case is somewhat intractable. Also, its realization in 
practice would probably require using data decisions, which might be 
unreliable when recovering from severe fades. Each of the above three 
strategies, by contrast, could be implemented using simple power 
measurements at the combiner output. The first two would involve a 
single, full-channel power measurement to drive the search over 6, 
and #2, while the third would involve three or more inband power 
spectrum density measurements, as described in Ref. 15. 

We emphasize that the third strategy is not optimal in the presence 


Aeq( fy= (7) 
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of equalization. Indeed, both of the first two strategies will be shown 
to yield superior outage statistics. Nonetheless, we will continue to 
refer to the gain pair (Ai, Bo) as optimal, but with quotation marks 
added. Thus, the gain pair (f, 82) is called optimal when there is no 
equalization and “optimal” when there is. 

Why even consider the third strategy in the presence of equalization 
if it does not yield the best detection? Because it serves another 
beneficial purpose, namely, sharply limiting the signal dispersion as 
seen at the input to the equalizer. The virtue in doing this is that it 
simplifies the circuitry needed to approximate the performance of 
ideal adaptive equalizers. For example, a digital equalizer might require 
fewer taps and/or lesser quantizing resolution by operating on input 
signals with less dispersion. Moreover, Gersho has shown that equal- 
izer convergence speed improves as input signal dispersion decreases.” 

Table I summarizes this section by giving, for each combination of 
space diversity and equalization to be considered, the assumed criteria 
for adapting (6;, 8). 


IV. ANALYSIS 
4.1 Receiver detection measure 


In current digital radio systems, the symbol rate (1/T) is typically 
75 percent of the channel bandwidth (W) and a is typically between 
0.25 and 0.5. Our simulations will assume typical symbol rates but, for 
convenience, will show results for a = 0. This case is especially easy 
to treat because the cosine roll-off function, C(/), reduces to a rectan- 
gle on [—1/2T, 1/2T] and this simplifies the analysis and computation 
of ISI power. At the same time, this approach leads to slightly 


Table |—Criteria for adapting (81, 82) under different space diversity/ 
equalization approaches 


Equalization 
Space Diversity None Reciprocal Equalizer 
None (B:, Be) a (1, 0) (8,1, Be) = (1, 0) 
Selection diversity (81, B2) = (1, 0) or (0, 1), (81, Bz) = (1, 0) or (0, 1), 
whichever yields larger pp whichever yields larger pp 
Combining diversity Optimal solution:* 1. B, = 1+ 0; | 8.| and 
(81, Bz) = (Bi, Be) Arg{8.} set to maximize 


combiner output s/n 

2. Bi = 1 +70; |B2| = 1; 
Arg{8.} set to maximize 
combiner output s/n 


3. “Optimal” golytion:* 
(81, B2) = (81, Be) 


* (81, Bo) is that gain pair for which pp is maximized when there is no equalization. 
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pessimistic results since, for given 1/T', ISI distortion decreases some- 
what with a.”?*! Our assumption for a, then, will serve both to simplify 
the analysis/computations and to yield worst-case estimates of detec- 
tion performance. 

To begin, we define 


Ps &4 The squared signal sample (excluding ISI and noise) in 
either baseband stream when a data value of +1 or —1 is 
being detected; 


P;, Py & The mean square ISI and noise, respectively, associated 
with the periodic samples in either baseband stream; 


and 
pp & Ps/(P1 + Py). (8) 


The latter is the detection signal-to-distortion ratio and can be related 
to an upper bound on the bit error rate via” 


Bit Error Rate < 2 exp(—pp/2). (9) 


This upper bound is particularly conservative when ISI dominates 
because this distortion component is generally peak limited, in con- 
trast to Gaussian noise. 

We first consider a receiver with no adaptive equalization (H.,(f) 
= 1+,0, in Fig. 1). For aw = 0, pp can be shown to be 








_ 3 
Pp >= M-1 
| Max(| |") (No ae 
[|H|* - Max(|H|*)] + (| Bi |? + | B2|?)/CNRJ Equalizer), 
where* 
H = { H(f)e*"dfT; (11) 
|H/? = + | H(A) |PofT: (12) 


CNR is the unfaded carrier-to-noise ratio per diversity branch in a 
bandwidth 1/T (typically close to 10°, or 60 dB); H(f) is given by (1); 
and ¢t, is the timing epoch. We assume the latter to be optimal, i.e., 
that value for which pp is maximized. This assumption is implicit in 
(10), where ¢, is specified to maximize | H |?. 


* In this analysis, all integrals have limits —1/2T and 1/2T. For convenience, we omit 
them from the equations. 
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As for H, it represents a complex gain for the signal vector sampled 
in each data interval. The imaginary part represents interference into 
each baseband rail from the cross-rail data in the same interval. 
Optimal carrier recovery amounts to multiplying H by a phase factor, 
e/° that makes it real and thus eliminates this interference. The 
resulting @ is the optimal carrier phase and the resulting signal gain is 
| H|, as used in (10). In the case of optimal combining diversity, the 
optimal @ is implicitly realized in the course of optimizing 8, and Bo. 
For both nondiversity and selection diversity, we will assume an 
optimal carrier phase without specifying how it is achieved. 

Now we consider a receiver with an ideal reciprocal equalizer. The 
signal-restoring property of this equalizer is such that only thermal 
noise enhancement modifies pp from its unfaded value. Thus, for the 
assumed system, | 


3 CNR (Reciprocal 


09 SS aaa i 
M-1 ; ; dfT Equalizer). 
(lal + 16) J oop 


The receiver detection measure we will use, and for which we will 
find probability distributions, is the bracketed quantity in (10) or (13). 


Thus, we define 
pa (Mo ) po (14) 


which can be viewed as a normalized signal-to-distortion ratio. It is 
a function solely of the response pair (H,(f), Ho(f)); the gain pair 
(81, B2); and the system parameters T' and CNR. Combining (14) with 
(9), we can upperbound the I required in an M-level system to achieve 
a specified bit error rate, BER,: 


(13) 





r, = 2(M — 1)In(2/BER,). (15) 


Decibel values of I, for various combinations of M and BER, are given 
in Table II. We will use these values later in assessing the simulation 
results. 


Table I|—Decibel values of I, for 
various M and BER, 
M 
BER, 16 64 256 
1073 18.81 25.04 31.11 
107* 19.96 26.19 32.26 
107° 20.87 27.10 33.17 
10-6 21.62 27.85 33.92 
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4.2 Formulas for 8, and 8 


We now present the appropriate gain pair (1, G2) as a function of 
(H,(f), H2(f)) for each of the various receiver processing approaches. 

In the cases of nondiversity and selection diversity, the specifica- 
tions are simple: For both equalized and nonequalized receivers, 


ij (1, 0); (No Diversity) 
(61, Be) — | (1, 0) or (0, 1), (16) 
whichever maximizes I’; (Selection Diversity). 


In the case of optimal combining diversity without equalization, 
the analysis is a bit complicated. The initial step is to rewrite (10) as 


_ 3 
PD M—1 


Ty \2 
. [Max Ms | ir (10’) 
ts L@r62) ([| 1]? — | H|7] + (1611? + [62|°)/CNR 

where |H|? and |H|* are quadratic in 6, and 2 and |H|? is, in 
addition, a function of t,. The maximization over 8, and $2 can be 
done analytically using the method outlined in Ref. 25. Briefly, the 
denominator (quadratic in 8, and $8.) is minimized subject to the 
constraint that H (linear in 6, and $2) equals 1 + j0. This problem 
can be solved using Lagrange multipliers and leads to the following 
result: 





B, = NUM,/A; — B2 = NUM,/A, (17) 

where 
NUM, = Hie’** — CNR(Hie/*" H¥H. — H¥e/ | H.|?); (18) 
NUM: = Hie/*" — CNR(Hie’" H,Hé — H¥e'*" |H,|?); (19) 
A = NUM, + NUM,; (20) 


and (_ ) is the average over f from —1/2T to 1/2T. (Because the 
processing is linear, an equivalent form of the solution is 8, = 1 + j0 
and Bo = NUM,/NUM,. For convenience only, we will use the solution 
as presented above.) 

These solutions for 8; and 82 maximize pp for a given value of t,. In 
our computations, the above formulas were used to obtain the inner 
maximum in (10’), and then a numerical search over t, was performed 
to find the outer maximum. 

Finally, we consider the use of combining diversity with adaptive 
equalization. As noted in Section 3.3, three distinct strategies for 
adapting 6, and #8 have been considered. Mathematical formulas for 
G, and 62 under these strategies are as follows: 
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1. 6; = 1 + j0 and both the magnitude and phase of 82 are adjusted 
to maximize the ratio of signal power to noise power at the combiner 
output. The solutions are 





Arg{62} = tan” [—Im(HH2)/Re(HiH2)] (21) 
and 
(62| = 2 ee et, (22) 
where 
Ya = | H2|? — | Hi |? (23) 
and 
Ys = [Re(HfH2)]’ + [Im(H¥H.). (24) 


2. 6B, = 1+ 70 and |6.| = 1, with the phase of 62 adjusted to 
maximize the same ratio. The solution for Arg{6.} is again (21). 

3. B, = B, and B. = Be, (17) through (20), where ¢, is that timing 
epoch for which pp in (10’) is globally maximized. (Again, 6; = 1 +j0 
and Bo = 82/8, would yield the same results.) 


4,3 Channel outage probabilities 


Our simulations lead to probability distributions for I over the 
ensemble of dual-channel response pairs. The probability distribution 
obtained for a given receiver technique will be denoted by P(T), 
representing the expected fraction of multipath fading seconds for 
which the receiver detection measure lies below IT. Thus, if I, is the 
value that I must exceed for acceptable performance, then P(T,) is 
the conditional probability of channel outage (i.e., conditioned on 
the occurrence of fading). Further, if P,.g is the maximum permissible 
value of that probability for meeting the system outage objective, the 
receiver technique can be considered adequate if and only if 
P(T.) S Preg- 

To apply this interpretation to our simulation results requires 
numerical data for both I, and P,.q. The former is derived in Section 
4.1 and quantified, for various combinations of M and BER,, in Table 
II. Derivations of P,.q are outlined in the Appendix and quantified in 
Table III. These results are for short- and long-haul systems in the 
4-,6-, and 11-GHz bands, both with and without protection switching. 


4.4 A performance measure for equalized receivers 


A second performance measure for which we shall obtain probabil- 


920 TECHNICAL JOURNAL, APRIL 1985 


Table II|—Values of Preq for various system conditions 


Protection Switching 
None One Protection 
Band System Length (Unprotected) Channel 

4 GHz (W = 20 MHz) Long haul 2.84 x 107 6.88 x 107° 
Short haul 4.54 x 107? 2.75 X 107? 

6 GHz (W = 30 MHz) Long haul 1.89 x 10-* 6.88 x 1073 
Short haul 3.03 x 107° 2.75 < 107? 

11 GHz (W = 40 MHz) Long haul 1.03 x 10~* 5.87 x 107° 
Short haul 1.65 x 107° 2.35 X 107? 


ity distributions is the range of |H(f)| over the Nyquist bandwidth, 
[—1/2T, 1/2T]. Thus,* 


Max! Hf) 
ip ea eaten < 1/27. 
oe {ahaa ik AS abo ee) 


The quantity R is of interest in cases where adaptive equalization is 
used, for reasons given in Section 3.3. It is this index of input signal 
dispersion that has been shown to provide a measure of convergence 
speed.*” We will see that its probability distribution, P(R), depends 
strongly on diversity approach. 


V. SIMULATION STUDY AND RESULTS 
5.1 Simulation approach 


A computer program was written that simulates response pairs 
(H,(f), He(f)) in accordance with the statistical model discussed in 
Section II and, for each simulation trial, analyzes the receiver tech- 
niques described in Section III using the methods of Section IV. 

The system parameters studied were channel bandwidth (W) and 
unfaded carrier-to-noise ratio (CNR). As discussed in Section 4.1, the 
symbol rate was taken to be 75 percent of the bandwidth (i.e., 1/T = 
15, 22.5, and 30 megabaud, respectively, for W = 20, 30, and 40 MHz) 
and the cosine roll-off factor was taken to be zero. The values used 
for CNR ranged from 51 to 67 dB, with 63 dB taken to be the most 
“typical” value. The channel parameters studied were dt, p and p, as 
discussed in Section 3.4, with the “typical” values taken to be 4 ns, 
0.5, and 0.65, respectively. 

For each combination of receiver technique and parameter set, either 
or both of the quantities [ and R were computed for each of four 


* We use the Nyquist bandwidth in defining this measure, noting that it is the 6-dB 
bandwidth of the system response for any value of the cosine roll-off factor. 
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thousand response pairs. From the population of quantities obtained 
for each such “simulation run,” we then computed the appropriate 
probability distribution(s), i.e., P(T) and/or P(R). 

The response pair for each trial was obtained via Monte Carlo 
generation of the a’s and b’s in (4), according to the assumed model. 
The same set of four thousand response pairs was used in each 
simulation run, except for minor differences related to changed values 
for 6t, p, and p. Our purpose was to permit comparisons among receiver 
techniques for the same fading channel ensemble. 

We confirmed, in two ways, that the ensemble of four thousand 
response pairs generated suitably represents the model, namely, (1) by 
computing the ensemble statistics of a, az, |b, |, |b2|, etc., and com- 
paring them with those predicted by the model; and (2) by computing 
P(YT) and P(R) for several different sets of response pairs and exam- 
ining the agreement between them. We thereby established that the 
simulations are accurate and consistent down to probability values of 
0.005 or lower, which is sufficient for our purposes. 


5.2 Results for combining diversity with no equalization 


This is the receiver technique of greatest interest to this study, as 
we wish to determine the multipath-combatting potentialities of those 
routes equipped with space diversity antennas. We computed P(I) for 
this case over a wide range of system and channel parameters. 


5.2.1 System parameters: W and CNR 


Figure 2 shows curves for P({) with W as a parameter. All other 
parameters are assumed to have “typical” values, i.e., CNR = 63 dB, 
6t = 4 ns, p = 0.5, and p = 0.65. The fairly strong influence of W is no 
surprise. The benefit gained by linearly combining H,(f) and H,.(f) 
depends on how effectively this combining cancels their dispersive 
components. Cancellation is most readily achieved when these com- 
ponents are linear in jw over the bandwidth, a condition most nearly 
approximated at smaller W.’° As W increases, higher-order terms in 
jw become prominent and a linear sum of two responses is less able to 
cancel the dispersions. In this light, we see that three or more diversity 
branches might yield significant gains at larger bandwidths. 

To apply the results of Fig. 2 to a concrete example, we assume a 
protected short-haul 64-QAM system with a BER objective of 107°. 
Table II shows that I’, ~ 28 dB for this case, and Table III shows that 
P(T,) must lie below 0.0275, 0.0275, and 0.0235 for W = 20, 30, and 
40 MHz, respectively. We thus see that using optimal combining 
diversity without equalization would suffice to meet outage objectives 
at 4, 6, and 11 GHz (the latter case being somewhat marginal). To 
meet long-haul objectives, however (see Table III), or to operate 
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CNR =63 dB: Ot =4 ns; p = 0.5; 9 = 0.65 
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Fig. 2—Probability distributions for I in receivers with optimal diversity combining 
and no equalization. The parameter is channel bandwidth, W. All other system and 
channel parameters have “typical” values: CNR = 63 dB, 6t = 4 ns, p = 0.5, and p = 
0.65. 


without protection switching or at 256 levels, optimal combining alone 
would be marginal or downright inadequate, at least for 6- and 11- 


GHz systems. 
To study the influence of CNR, Fig. 3 shows P(T) for W = 30 MHz 
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CNR = 51 dB 


P(T) 


NO EQUALIZATION, OPTIMAL COMBINING 
W = 30 MHz; Ot =4 ns; p =0.5; p = 0.65 





20 25 30 35 40 
T IN DECIBELS 


Fig. 3—Probability distributions for I in receivers with optimal diversity combining 
and no equalization. The parameter is CNR. Channel bandwidth W = 30 MHz, and all 
other system and channel parameters have “typical” values. 


with CNR parameterized over the range from 51 to 67 dB. The trend 
is similar for all other W. The sensitivity of the results to CNR is seen 
to increase as this quantity decreases; this is because noise becomes 
an increasingly significant factor (as opposed to intersymbol interfer- 
ence) as CNR decreases. 
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Fig. 4—Probability distributions for I in receivers with optimal diversity combining 
and no equalization. The primary parameter is 6t, with results given for W = 20 and 40 
MHz. All other system and channel parameters have “typical” values. 


5.2.2 Channel parameters: dt, p, and p 


Figure 4 shows P(T) for 6t = 2, 4 and 6 ns for each of two bandwidths 
(W = 20 and 40 MHz). Not surprisingly, the sensitivity of performance 
to this parameter is greater at larger bandwidths. In the remainder of _ 
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Fig. 5—Probability distributions for I’ in receivers with optimal diversity combining 
and no equalization. The vertical bars span the range from minimum to maximum over 
the p-range [0, 1.0] and the solid curve is for p = 0.5. Channel bandwidth W = 30 MHz 
and all other system and channel parameters have “typical” values. 


this study, we will confine ourselves to the assumed “typical” value of 
4 ns. 

Figure 5 shows the effect of p, the minimum-phase probability, when 
W = 30 MHz. The vertical bars show the range of P(IT) as p goes from 
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0 to 1.0, and the solid curve is for the assumed “typical” value of 0.5. 
Similar results apply for other values of W. The highest point on each 
bar is for either p = 0 (all Hi(f) and H.(f) minimum phase) or p = 
1.0 (all H,(f) and H2(f) nonminimum phase), and the lowest point is 
for p = 0.5. The explanation for these trends is simple: The conditions 
p = 0 and p = 1.0 provide, on the average, for the greatest shape 
similarities between H,(f) and H2(f), while the condition p = 0.5 
provides, on the average, for the greatest dissimilarities. In any event, 
Fig. 5 shows that a precise quantification of this parameter is not 
important for the situation under study. 

We also examined the effect of the correlation factor p, for W = 30 
MHz, over the range from 0 to 1.0 (we did not seriously consider 
anticorrelation). The spread was found to be very small, suggesting 
that the amplitude factors a; and az in (4) are of minor import when 
CNR is high and that, therefore, the channel parameter p need not be 
precisely known. Clearly, then, poor performance occurs primarily 
when the relative variations of H,(f/) and H.(f) are at once highly 
dispersive and (by random chance) highly similar. We invoked this 
observation to derive P,., in the Appendix. 


5.3 Comparisons with non- and selection diversity 


Figure 6 shows P(T) for non-, selection, and optimal combining 
diversity without equalization. These comparisons are given for two 
bandwidths (W = 20 and 40 MHz). The parameters CNR, 6¢, p, and p 
are assigned their “typical” values. For non- and selection diversity, 
the values of 5t and p have no influence on P(T). The values of CNR 
and p were found to have some impact (the statistics improve slightly 
with increasing CNR for both non- and selection diversity, and im- 
prove slightly with decreasing p for selection diversity), but we did not 
study these relationships in detail. 

For both bandwidths, receiver performance without space diversity 
or equalization is seen to be disastrous, as expected. Selection diver- 
sity improves matters noticeably, but is still far from adequate. The 
improvement in going to optimal combining diversity is seen to be 
dramatic. 


5.4 Results for receivers with equalization 


For receivers using ideal reciprocal equalizers, we have obtained 
both P(T) and P(R). We can thus examine trade-offs between detec- 
tion performance (measured by P(T)) and dispersion at the equalizer 
input (measured by P(R)). For all diversity arrangements, both P(T) 
and P(R) depend weakly (if at all) on 6¢ and p. In addition, P(T) 
depends very weakly on W and P(R) depends very weakly on CNR. 
In several cases, the correlation factor p was found to have a discernible 
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Fig. 6—Probability distributions for I in receivers with no equalization. Results are 
for non-, selection, and optimal combining diversity, for W = 20 and 40 MHz. All other 
system and channel parameters have “typical” values. 


impact. This is because ideal equalizers tend to capitalize on total 
signal power and lower values of p (greater independence between a, 
and a2) correspond to statistically higher values of total signal power. 
For this part of the study, we confined ourselves to the “typical” and 
somewhat conservative value of 0.65. 
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Fig. 7—Probability distributions for I in receivers with ideal reciprocal equalization. 
Results are for non-, selection, and combining diversity, with three gain-setting algo- 
rithms considered for the latter. Channel bandwidth W = 30 MHz and all other system 
and channel parameters have “typical” values. 


5.4.1 Probability distributions for T 


Figure 7 shows P(I) for non-, selection and combining diversity, 
with three distinct gain-setting algorithms considered for the combin- 
ing case. Some observations are as follows: 
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1. While the results shown are for W = 30 MHz, they would barely 
change for other bandwidths. With this in mind, we can see from Figs. 
7 and 2 that a reciprocal equalizer alone would be less effective than 
an optimal diversity combiner alone for W = 20 MHz, but more 
effective for W = 30 MHz. 

2. For reciprocal equalizers, I’ is related strictly to thermal noise 
enhancement. Thus, this quantity scales directly with CNR (63 dB in 
Fig. 7), that is, for every decibel of reduction (or increase) in CNR, all 
curves would shift by one decibel to the right (or left). 

3. For combining diversity, the best results shown are for the 
(61, B2) pairs that maximize combiner output power. Controlling both 
amplitude and phase is only slightly superior to controlling phase 
alone, which is only slightly superior to using selection diversity. Using 
the gain pair (f,, 82), while optimal in the absence of equalization, is 
seen to be fourth best in this case. We now examine the potential 
benefit of using these suboptimal gains. 


5.4.2 Probability distributions for R 


Figure 8 shows P(R) for W = 20 and 40 MHz. For each bandwidth, 
the results shown are for nondiversity, selection diversity, maximum- 
power combining diversity (phase-only adjustments), and “optimal” 
combining diversity, (81, 82) = (81, B2). The disparities among ap- 
proaches are seen to be strong, the best results corresponding to 
“optimal” combining. 


VI. CONCLUSION 


The results obtained in this study are subject to the usual uncer- 
tainties associated with finite simulations and a less-than-universal 
statistical model of dual-channel fading. Nevertheless, some general 
conclusions are possible, as follows: 

1. The use of appropriate space diversity combining (e.g., the scheme 
described in Ref. 15) could eliminate the need in some links for post- 
combiner equalization. 

2. In links where adaptive equalization 1s used, such a scheme could 
sharply reduce the signal dispersion at the equalizer input, thereby 
simplifying the equalizer design. 

3. For purposes of assessing various receiver approaches, the dual- 
channel multipath fading model reported by Rummler is, for the most 
part, well specified. The precise value of the minimum-phase proba- 
bility (p) appears to be of minor consequence, as is the precise value 
of the correlation factor (p) between the decibel values of a; and ag. 
The delay parameter (6t) introduced here could have a modest impact 
if larger than 4 ns, particularly for the larger channel bandwidths. 

4, Finally, the numerical results affirm that space diversity combin- 
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Fig. 8—Probability distributions for R for several diversity approaches, for W = 20 
and 40 MHz. All other system and channel parameters have “typical” values. 


ing and adaptive equalization in tandem comprise a formidable com- 
bination of receiver techniques in microwave digital radio. 
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APPENDIX 
Derivation of Preq 


We consider a specific W-Hz radio channel carrying digital traffic. 
It may be protected or unprotected and it may have dual space diversity 
or no diversity. Discernible multipath fading is assumed to exist, on 
one or both of the possible diversity branches, for T, seconds per hop 
per year, with T, given by (8). We define the following probabilities 
within that subset of time: 


P, & The probability of outage per hop as experienced by a user 
of the channel; 


P, 2 The probability of outage per hop as experienced by the 
channel itself (note that P, < P,; when protection switching 
is used); 


Preq & The value of P,; required to achieve a specified value of P,. 


We now make the following assumptions: 

1. The yearly outage objective on a total system route is 0.005 
percent (one way), apportioned uniformly among the route hops. 

2. The path length and terrain factor for each hop are taken to have 
average values: D = 25 miles and c = 1 in (3). 

3. Within any protection switching section, discernible multipath 
fading occurs on at most one hop at a time. 

4, The subset of multipath fading seconds on each hop is congruent 
over all W-Hz channels within the same common carrier band. 

5. Within that subset of multipath fading seconds, outage events 
are statistically independent from one W-Hz channel to another. 

The third, fourth, and fifth assumptions are germane to calculations 


SPACE DIVERSITY AND EQUALIZATION 933 


of Preq for systems with protection switching. The third assumption is 
slightly liberal, in that it leads to overestimation of P,.q; the fourth 
assumption, by the same token, is somewhat conservative; and the 
fifth assumption is discussed later. 

In systems without protection switching, P, and P, are the same, 
and so Pye, is just the value of P, dictated by the first two assumptions. 
For a long-haul system (maximum length = 4000 miles), the outage 
budget will be met if each hop averages 10 seconds of outage per year; 
for a short-haul system (maximum length = 250 miles), that number 
is 160 seconds. Dividing each number by T, yields the required Preg. 
The results are shown in Table II for the 4-, 6- and 11-GHz bands. 

For systems with protection switching, the issue is more compli- 
cated, involving as it does the joint outage probabilities for channels 
within the same band. For most receiver techniques, outages in neigh- 
boring channels during multipath fading are not independent events, 
in which cases the fifth assumption cited above is too liberal. It might 
be quite valid, however, for the special case of optimal combining 
diversity without adaptive equalization. The reasoning is as follows: 
In such receivers, outage occurs primarily when H,(f) and H.(f) have 
similar shapes and, thus, dispersion cannot be reduced without exces- 
sive signal loss. There is no reason, however, why such a similarity 
between two uncorrelated responses would occur in two frequency 
channels at the same time, except by chance; hence, the independence 
assumption. | 

The case of optimal combining diversity without adaptive equali- 
zation is of major interest here, as we mean to explore the limits of its 
applicability in digital radio. Thus, we use the independence assump- 
tion in deriving P,.,, speculating that the result is accurate for that 
type of receiver and somewhat elevated for most others. 

To proceed, let N be the number of channels in a given band and 
assume that one of them is used for protection. Given an outage in a 
traffic-bearing channel, the probability that it will not find the pro- 
tection channel available is 


p= P+ (2 





P, + smaller terms. 


The first term is the probability that the protection channel itself is 
out; the second term is the probability that the protection channel is 
not out but one of the other N-2 channels is and switches first; and 
the smaller terms (neglected here) have to do with three or more 
channels being out at the same time. Outage for a given user occurs if 
its original channel is out (probability P,) and the protection channel 
is not available (probability P,). Thus, 
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N 
Py = PyP2 => Pi. 


We can now find the P, required to achieve a specified P,. The latter 
values, as before, are 10/7, and 160/T,, respectively, for long- and 
short-haul systems. The values of N are 12, 8, and 6, respectively, for 
the 4-, 6- and 11-GHz bands. The resulting values for P,., for each of 
various systems are shown in Table III. 
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A Controlled Impedance Robot Gripper 
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The design of robot grippers has traditionally involved pneumatics or gear 
reduction electric drives that do not allow for the control of the stiffness of 
the mechanism. The manipulation of fragile objects such as eggs or light bulbs 
requires the gripper to be able to close on the object with minimal impact 
forces, and yet maintain a static grip force sufficient to firmly handle the 
object. This paper describes a two-fingered controlled impedance gripper where 
the impedance is imparted by electrical means. The analog control system 
allows independent control of the effective mechanical mass and damping of 
each finger, as well as additional control of common-mode versus differential- 
mode response. The gripper is modeled by a computer simulation that consists 
of a set of nonlinear differential equations with time-varying feedback param- 
eters. In the final form the model has ten degrees of mechanical freedom and 
eight electrical poles, i.e., an 18th-order nonlinear differential equation. The 
dynamic equations are described and results given are compared with actual 
results obtained from laboratory experiments with the gripper. An Appendix 
describes the method of solution used for the nonlinear differential equations. 
The accuracy of the simulation has been verified by measurements on the 
actual gripper. This analysis has resulted in the design of a gripper control 
system capable of providing controlled compliance and reduced finger impact 
forces while maintaining a quick response and firm grasp. 


I. INTRODUCTION 


Many compliant mechanisms have been built using mechanical 
means for obtaining the compliant properties. For example, the Re- 
mote Center Compliant manipulator (RCC), developed by Whitney 
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and Nevins’ at Draper Labs, utilizes specially designed springs in a 
configuration that allows free angular movement about a radial center 
while maintaining high stiffness along the radius. This device is 
particularly useful for inserting a pin in a hole or for analogous tasks, 
but the amount of compliance is fixed while in operation. 

Similar devices have been developed by the Lord Corporation and 
others. One interesting device was designed by M. Cutkowski.” It 
replaces the fixed springs of the RCC device with fluid-filled elastomer 
bubbles. The compliance thus formed is nonlinear and adjustable by 
changing the fluid pressure. This feature makes the compliant wrist a 
potentially active device with bandwidth limited only by the rapidity 
with which fluid pressure may be changed. 

Another compliant manipulator was designed by Hanafusa and 
Asada, who developed a three-fingered gripper with spring-driven 
fingers.° The fingers are arranged in a circle and have rollers on the 
fingertips. Much of their work has involved determining stable grip- 
ping orientations for irregularly shaped objects. 

Dynamic computer control of the amount of compliance and damp- 
ing imparted to the gripper is desirable since this allows the gripper 
to adapt to various grasping tasks as they proceed without the need to 
exchange mechanical parts such as springs.* For a two-fingered grip- 
per, a “common-mode” compliance and “differential-mode” compli- 
ance can be made usefully different so that a gripper may use a delicate 
grasp and still hold an object firmly in place. : 

Figure 1 illustrates the construction of the two-fingered gripper. 
The fingers are mounted on linear bearing slides and driven through 
a rack and pinion by two dc servomotors. The two fingers are driven 
independently by two dc servomotors. Thus, the gripper can close 
from either side and need not be centered on an object before closing. 

A particularly important feature of this mechanical design is the 
rack and pinion servomotor device. This mechanism is mechanically 
back drivable, meaning that a force applied to the finger will cause the 
motor to move. A typical worm drive mechanism, for example, would 
not have this characteristic because of high internal friction. This 
feature allows the motor drive mechanism to function as a mechanical 
impedance (under computer control) applying all of the appropriate 
effects of mass, damping, and compliance directly to the fingers. 

To control the position of the gripper fingers, a means of measuring 
their positions is necessary. A capacitive measuring system was devel- 
oped (see Ref. 5) that has the properties of low noise 
(s/n > 90 dB), low mass, and fast transient response. It consists of a 
metal rod sliding inside an insulated metal tube. The capacitance of 
the rod and tube varies as a function of the position of the finger. A 
capacitance measurement circuit produces a voltage proportional to 
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Fig. 1—Instrumented controlled impedance gripper. 


the capacitance. The position measurement error developed by this 
system is less than 1 percent, while the step response rise time is 
about 150 us. 

The dc servomotor performance can be modeled by an electrical and 
mechanical system consisting of an armature resistance and induct- 
ance, a back Electromotive Force (EMF), armature inertia, and bearing 
friction. This forms a third-order differential equation describing the 
current-voltage relationship as well as the torque output of the motor. 
Since torque is a quantity that must be controlled and torque is an 
essentially linear function of motor current, it is desirable to drive the 
motor with a current source. This reduces the motor equations to 
second order and greatly simplifies the design of a controller. Motor 
voltage is then described by an auxiliary equation. 

All dynamic systems can be represented by a set of states or state 
variables which describe the system energies.® The energies can be 
divided into two basic classes: static/potential and dynamic/kinetic, 
or analogous terms. For example, the motor state can be described by 
armature position and rotational velocity. State equations are typically 
written in matrix notation. In particular, most forms of dynamic 
systems can be represented by two matrix equations: 


x(t) = A(x, ¢) x(t) + B(x, ¢) u(t) (1) 
y(t) = C(x, t) x(t) + D(x, ¢) u(Z) (2) 
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where x is the state variable vector and (1) is the “state equation.” 
Matrix A(x, ¢) describes the feedback part of the system and matrix 
B(x, t) describes the feed forward part from the input vector u(t). 
The “output equation,” (2), describes all auxiliary variables that may 
depend on the system state or input vector but do not affect the system 
state. 

As an example, let us model the motor described earlier in state 
variable form, with a voltage source as the input. The mechanical state 
variables are motor shaft position and velocity, and electrical energy 
is best described by the current flowing in the inductance of the 
armature. Given these three state variables, the coefficients of the 
motor equation are 


F 1 0 | 
A=|0 -P/M K,/M (3) 
0 —-K./L —R/L 
0 
B .| 0 | (4) 
1/L 


where K; is the motor torque constant, M is the armature moment of 
inertia, P is the viscous friction coefficient, K. is the inducted back 
EMF coefficient (EMF = K, X Vel), L is armature inductance, the 
state variable vector is 


armature position (6) 
x = (5) 


armature velocity (6) 
armature current I 


and the input u(t) = V,(t). The motor torque T is a direct function 
of the armature current, giving for the output equation 


T = y = Cx(t) + Du(t) = [0 0 K,)]x(é), (6) 


where D = 0 in this case. 
The third equation in (8) can be eliminated by driving the motor 
with a current source. The coefficients of the state equation become 


0 1 
oul : — 7) 


0 
B= a (8) 


where the state variable vector consists of armature position and 
velocity only, and the input u(t) = J,(t). The output equation becomes 
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T = Du(t) = Kil;(t), (9) 


where C = 0 in this case. 

The control system is described in the next section. Details of the 
computer simulation are discussed in Section III. In Section IV the 
results of various parameter studies are presented. Section V describes 
some experiments with the completed gripper system. Finally, some 
concluding remarks concerning certain significant results are offered 
along with some plans for future study. 


Il. GRIPPER CONTROL SYSTEM 


The principal function of the controller is to position each gripper 
finger appropriately for approaching an object and to close the gripper 
fingers on the object. Another goal of the gripper control system design 
is to provide a means for affecting the mechanical properties so that 
gripper compliance and dynamic performance can be modified under 
computer control. For example, it is desirable to minimize impact 
forces due to the inertia of the gripper finger and motor armature 
masses while allowing rapid response in the mechanism. One way to 
reduce impact forces is to reduce the mass of the fingers. Since the 
mechanical mass of the mechanism has already been minimized by 
design, it is desirable to further reduce the effective mass of the system 
by appropriate control variable feedback. 

The effective mass of the gripper mechanism can be controlled 
electrically by feeding back an acceleration term in the control system. 
It is well known that effective damping in the system is controlled by 
velocity feedback. This was the basis for an early design of the gripper 
control system. Subsequently, after nonlinear computer analysis of 
several system configurations, a system design evolved which is de- 
picted in Fig. 2. 

The control system for the gripper consists of an analog section and 
a digital interface. Analog feedback methods were chosen after nu- 
merous computer simulations indicated that a microprocessor would 
not be fast enough to perform the complicated control calculations 
necessary at a sufficiently high sample rate. The control circuitry 
consists of 24 operational amplifiers, six 8-bit Multiplying Digital-to- 
Analogue Converters (MDACs), two 10-bit DACs, and a digital inter- 
face. The digital interface communicates with a SUN 68000 microcom- 
puter running under New Real-Time Executive (NRTX) operating 
system.’ Mechanical compliance is controlled by position feedback, 
damping by velocity feedback, and effective mass by acceleration 
feedback. The velocity and acceleration signals are obtained from the 
position signal by electrical differentiation. In order for this technique 
to work well, the position signal must have very low noise. The 
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Fig. 2—Gripper control system. 
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amounts of position and velocity feedback are controlled by computer 
through the digital interface and 8-bit MDACs indicated by coeffi- 
cients b and c in Fig. 2. Finger positioning stiffness is controllable 
from zero to about 23,000 N/m in increments of about 90 N/m. 
Acceleration feedback gain is fixed by a trimmer potentiometer and is 
set as high as possible consistent with stable operation of the gripper. 
Command positions [r;(t) and re(t)] are controlled by 10-bit DACs 
for a positioning resolution of about 40 ym. 

An additional feature of the control system is a means of cross- 
coupling the two finger positions so that an effect on common-mode 
compliance is obtained. Common-mode damping is provided by com- 
mon-mode velocity feedback. Both common-mode position and veloc- 
ity gains are controlled by 8-bit MDACs (coefficients e and f in 
Fig. 2). This common-mode feedback allows the gripper fingers to, for 
example, hold an object with a soft grip (differential mode) while 
stiffly holding the object’s position with respect to the gripper body 
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(common mode). Alternatively, the common-mode stiffness can be 
made zero so that the finger positions can conform to an off-center 
object when they are closing (such as a peg in a hole). 

Common-mode compliance control can be illustrated by a few ex- 
amples. Ignoring velocity and acceleration gains, the feedback equa- 
tions for the gripper fingers are 


I, = ay(ry — x1) + f (x1 + x2) (10) 
In = ao(Te — X2) + f(x1 + x2), (11) 


where /, is the motor current control signal for motor 1 (which may 
be scaled by fixed gain constants), a; is position gain associated with 
finger 1, f is the common-mode position gain, r, is the commanded 
position for finger 1, x; is the actual position of finger 1 (positive 
direction is motion toward the center of the gripper), and x2 is the 
actual position of finger 2, etc. The gripper fingers always tend to 
move so that J; and J, are zero. The amount of deviation of the values 
of J, and J, from zero is proportional to the reactive force generated 
by the servomotor. 

By appropriate choice of a;, az, and f, several different modes of 
operation can be obtained. Under normal operating circumstances, the 
values for finger 1 (a, b, and c) are chosen to be the same as those of 
finger 2 for balanced operation. If f is set to one half the value of a, 
for example, then the common-mode stiffness is very high and differ- 
ential-mode stiffness is zero. The control equations become 


I, = ar, + a(xe — X1)/2 (12) 
Ig = arg + a(x1 — X2)/2, (13) 


where a = a, = a2. This results in an effect where each finger 
individually has no stiffness with respect to the gripper body as long 
as the other finger is free to move. However, if both fingers are forced 
to move in the same direction, considerable stiffness is encountered. 
In this mode the fingers tend to stay equally distant from a position 
value centered between the two-finger command values r; and ro. If 
the values of the a coefficients are increased, a soft grip can be 
established with very firm centering of the object. This is useful when 
the robot arm is in motion as it keeps the fingers and object centered 
without unnecessary gripping force. 

Alternatively, if the a coefficients are made zero and f is negative, 
common-mode stiffness is nullified and all stiffness appears as differ- 
ential-mode stiffness. The control equations reduce to 


I, = —a(x, + Xa) (14) 
I, = —a(x; + Xo). (15) 
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This means that the gripper fingers are stiff with respect to each 
other, but exhibit no stiffness with respect to the gripper body. This 
mode of operation is particularly useful in picking up objects that may 
not be positioned exactly in the center of the gripper. The fingers will 
easily conform to the position of the object without applying significant 
lateral forces. Varying amounts of common-mode versus differential- 
mode stiffness can be obtained with values of f between zero and one- 
half a. 

We will introduce the state feedback control equations by first 
considering a single finger with position, velocity, and acceleration 
feedback with a current-driven servomotor. Initially, we will assume 
that the velocity and acceleration are available directly from the finger 
by appropriate sensors (later we will remove this requirement). We 
will also assume a linear system until later when nonlinearities are 
added. Then we can write a feedback equation for the motor current 
as follows: 


I,(t) = —ay(t) — by(t) — c[y(t) — r(t)], (16) 


where r(t) is the reference position input function and a, b, and c are 
gain coefficients. Substituting (16) into (7) through (8) and combining 
like terms yields 


0 0 
oa =. + aK,) —(P + bK,)/(M + ” 0 


0 
ca Lace! + ak on 


where, again the state vector, x, consists of finger position y and 
velocity y, and the input function u(t) = r(¢). 

The lower left element of A is now nonzero, indicating that a 
position-sensitive term is present. This term is effectively a spring 
constant divided by a mass. For comparison, examine the state equa- 
tions for a classical elementary mass, spring, and damper problem: 


0 1 
A = _ as * (19) 


0 
B= F Ha (20) 


where k is the spring constant, p is the damping coefficient, and m is 
the mass. The input function u(t) applies a force to the mass, which 
is attached to a stationary support by the spring and damper. Note 
that in (17) through (18) the spring constant is replaced by (cK,), the 
damping by (P + bK;), and the mass by (M + aK;,). Thus, these 
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mechanical parameters can be independently controlled by the elec- 
trical feedback gains (a, b, and c) in the control system. 

The second finger is controlled in a similar manner. The two fingers 
thus operate independently in this case. Figure 2 also shows control 
signal passing from each finger to the other. These signal paths provide 
additional common-mode position and velocity feedback, which in- 
creases the common-mode stiffness. It was found, during computer 
simulation, that common-mode acceleration was not useful and that 
it tended to destabilize the control system. The complete feedback 
equation for motor drive current is 


[,(t) = -(a + d) y(t) — dyo(t) — (b + e)yi(t) — exe(t) 
—(c + f)yi(t) — fye(t) + (e + f)ur(t) + fue(t). (21) 


Closed-form solutions of the above equations are not readily obtain- 
able for a nonlinear system. The construction of our gripper involves 
servomotors with power limitations. The motor driver has both voltage 
and current limitations (which must be treated separately). The rack 
and pinion have backlash, and the mechanical parameters change 
abruptly when an object is grasped. The rack and pinion gear backlash 
has turned out to be a small effect, but it was not initially clear if this 
hysteretic effect would cause system instability with the electrical 
feedback. 


Hl. COMPUTER SIMULATION 


A fourth-order Runge-Kutta algorithm was used to study the time 
response of the nonlinear system just described. The algorithm used 
is modified somewhat from standard routines found in the literature 
and is described in Appendix A. The state equations must be expanded 
to include the various nonlinearities. After doing so, the two-fingered 
model requires an eight-dimensional state vector. Figure 3 is a diagram 
of the model and includes an object to grasp (in this case a rubber 
ball). For the time being, we will still assume that the finger velocity 
and acceleration are directly available. Later, we will study the influ- 
ence of electrical time constants on the mechanical performance. 

The model shown in Fig. 3 consists of five mechanically independent 
members. They are the two fingers, two motor armatures, and the 
rubber ball. Each mechanical member will require two state variables 
to describe its position and velocity. The ten-dimensional state vector 
is 

1. Left finger position 
2. Left finger velocity 
3. Left armature position 
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. Left armature velocity 

. Right finger position 

. Right finger velocity 

. Right armature position 
. Right armature velocity 
. Ball position 

10. Ball velocity. 

The state equation describes forces that exist on each mechanical 
member. These forces may be due to acceleration, frictional drag, or 
spring forces. Also included are forces applied to each member by 
another member when they are in contact. To introduce the state 
feedback matrix, let us first consider the case where no contact is 
being made between any of the mechanical members and the common- 
mode feedback is ignored. Then the coefficients of the state equation 
are 


Oo MAND oO > 


0 1 Oo 0 0 0 O 0 00 
0* —P/M* 0* — 0# 0 0 O 0 O* OF 
0 0 0 1 0 0 O 0 00 
—c K,/m* —b K,/m* 0* —p/m* 0 0 0 0 0 0 
0 0 0 0 0 1 O 0 00 
ies 0 0 oO 0 o* —P/M* 0*  0* O* 0* (22) 
0 0 0 0 0 0 0 1 0 0 
0 0 0 0  —c K,/m* —b K,/m* 0* —p/m*0 0 
0 0 0 0 0 0 O 0 01 
0* o* 0 0 o* o* 0 0 O* O* 
0 0 0 0 
0 0 0 0 
0 0 0 0 
K, 0 —a K,/m 0 
0 0 0 0 
a 0 0 0 0 (23) 
0 0 0 0 
0 K, 0 —a K./m 
0 0 0 0 
0 0 0 0 


where P is the damping coefficient of the finger slide mechanism, M 
is the mass of the finger, p is the damping coefficient of the motor 
armature, and m is the mass of the armature. The rubber ball rolls on 
a table with negligible friction. | 

The input function u(t) is a vector of four dimensions. The first 
two dimensions are the left and right finger position reference inputs. 
The last two dimensions are inputs from the implied left and right 
accelerometer sensors. This variation of the feedback equation is useful 
since acceleration is not one of the state variables available. Alterna- 
tively, the time derivative operator could be used in the A(x, t) matrix. 

The A(x, t) matrix becomes more complicated when any of the 
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Fig. 3—Model of two-fingered gripper. 


mechanical elements come into contact. Some or all of the elements 
of the matrix indicated with an asterisk in (22) become nonzero, and 
additional forces are coupled to the various mechanical parts in a 
symmetric manner. Equal and opposite forces are introduced into two 
rows of A(x, ¢) when contacts are made. Thus, the matrix form of 
(22) is helpful in illustrating the symmetry in the problem. For 
example, when the left finger contacts the ball, elements 1, 2, 9, and 
10 of row 4 are affected. At the same time, the corresponding elements 
of row 10 are also affected with identical magnitude, but with opposite 
signs. The motor drive nonlinearities are saturation limits on both 
voltage and current. These enter directly into the motor torque equa- 
tions. | 

The actual nonlinear state equation would consume more than a 
single page. A C program listing of the complete state equation is given 
in Appendix B. This listing of the state equation also includes the 
electrical time constants, which deserve some discussion. 
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Heretofore, we had implied that the acceleration of the left and right 
fingers would be available via some sensor, such as an accelerometer. 
In the actual gripper design, the acceleration and velocity are obtained 
by electrical differentiation of signal from the position sensor. The 
design of an electrical differentiator purposely incorporates a pole in 
the transfer function to limit the high frequency gain and thus reduce 
the amount of noise developed by or passed through the circuit. It is 
desirable to make this frequency limit as low as possible to minimize 
the circuit noise. However, the accuracy of the differentiation may 
suffer if the pole location is too low in frequency. An analysis of the 
effects of differentiator pole position was therefore necessary. The 
motor drive amplifier also has a dominant pole at a moderately low 
frequency (about 1500 Hz) and its effect is also included in the state 
equation. 


IV. RESULTS OF THE SIMULATION 


Transient response analysis was performed for various parameter 
settings to determine the effects of different amounts of feedback, the 
effects of various ball masses and compliances, and the effects of the 
differentiator time constants. A fourth-order Runge-Kutta algorithm 
was employed with modifications of the step size correction algorithm 
as described in Appendix A. The time response output is presented in 
graphical form that is easily read and gives a good indication of the 
performance. 

Several parameters were initially established for the analysis. These 
include the mechanical specifications for the dc servomotors and the 
finger mechanisms. The servomotors were selected for this application 
based on their physical size, torque capability, and the fact that they 
are completely characterized by the manufacturer. The published 
specifications for the motors are 


Armature resistance (Ra) 6.12 0 
Armature inductance (La) 0.00275 H 

Motor EMF constant (Ka) 0.0446 V-s/rad 
Motor torque constant (Km) 0.0446 N-m/A 
Armature inertia (m) 1.54 x 107° N-m-s?/rad 
Armature damping (p) 2.74107 N-m-s/rad. 


The symbols in parentheses refer to variables in the computer program 
(Appendix B). In addition, specifications for the rack and pinion gear 
mechanism used are 


Pinion gear radius (rg) 0.0081 m 
Pinion gear backlash (H) 508x107 m 
Finger mass (M) 0.6 kg 
Finger damping (P) 0.5 N-s/m. 
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The following parameters are variable in the simulation and can be 
set to establish interesting initial conditions: 
(a) Mechanical parameters: 


Initial ball position = (y[8]) 


Ball radius (R) 
Ball stiffness (k3) 
Ball mass (m3) 
Ball damping (p3) 


(b) Electrical parameters: 


RC time constants (RC) 
Position feedback gain for each finger (c) 
Velocity feedback gain (b) 
Acceleration feedback gain (a) 
Common-mode position gain (f) 
Common-mode velocity gain (e) 
Common-mode acceleration gain (d) 
Motor driver voltage limit (Vmax) 
Motor driver current limit (Imax). 


Several of these items introduce nonlinearities into the state equation. 
Two voltage programmable power supplies were used for motor drivers. 
These power supplies are linear units capable of operating either as 
voltage or current sources. The supplies used in this study are limited 
to 50 volts and/or 8 amperes. They can operate as either a power 
source or sink. They can also be controlled by either a voltage input 
or via a digital port, making them useful for microprocessor control. 

Initially, the transient analysis is conducted with step position 
command inputs. Later, we will discuss exponential input functions 
and their effect on performance. To perform a typical simulation the 
initial conditions are set so that the gripper is in the fully open 
position. Thus, the gripper is opened to a width of 6.35 cm. The step 
input functions are then applied and the gripper closes on an object, 
in this case the ball. Time responses of each finger and the ball 
position are computed along with the associated motor voltages and 
currents. Also calculated are the ball position (and velocity), the 
impact force being applied to the ball by each finger, and the estimation 
errors of the electrical differentiators. This information is all displayed 
on a single page in graphical form for each test response. 

The gripper was designed to allow each finger to move 3.8 cm. The 
fingers have about 1.25 cm of travel overlap, so in the fully opened 
position they are about 6.35 cm apart. The positional reference for the 
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Fig. 4—(a) Step response with position feedback. (b) Step response with position and 
velocity feedback. 


fingers was established so that the closed position of the gripper with 
the fingers centered is considered the zero position of each finger. The 
left finger opens in the minus x direction and the right finger opens 
in the plus x direction. 

Figure 4 illustrates the effects of power supply limits, gear backlash, 
and velocity feedback on the time response. For these responses the 
gear backlash was set to 1 mm to make the hysteresis effect obvious. 
This is roughly equivalent to removing every other tooth from the 
gear. Figure 4a shows a response with position feedback only. Without 
velocity feedback the mechanical damping is obviously inadequate. 
The gear backlash is only slightly visible in the plot of position versus 
time and appears as a slight jaggedness in the response. The voltage 
and current plots show the effects of the motor driver limits and the 
inductance of the motor. Note that initially the voltage is saturated at 
—50 volts. The current goes to nearly —8 amps but then decreases 
until the driver switches to +8 amps. This ramp effect is due to the 
induced back EMF of the motor as it gains velocity. The voltage 
response shows similar characteristics during current saturation. 

Figure 4b shows a response with viscous damping applied by velocity 
feedback. In this case the velocity gain is set to properly complement 
the position gain so that settling time of the finger is somewhat less 
than “critically damped.” This amount of damping results in minimum 
settling time to the final position since a very slight amount of 
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overshoot in the response (a few percent) results in much faster 
convergence to the desired control point. The effects of gear backlash 
can now be seen in the voltage and current plots. During the dynamic 
part of the response there is very little evidence of gear backlash. The 
gear is first in contact on the forward thrusting (accelerating) face and 
then on the decelerating face. After the position of the finger becomes 
stable, the gear must continually switch faces to maintain the finger 
in position. This backlash effect shows up as a small oscillation in the 
drive current and voltage. The actual gear backlash of the gripper is 
much smaller than this and its effect will not be apparent in the time 
responses that follow. 

When the fingers come into contact with the ball, impact energy 
and forces are exerted that may cause damage to the ball (or object) 
and/or gripper mechanism. These impact forces can be affected sig- 
nificantly by acceleration feedback. Figure 5 illustrates this effect. A 
typical damped closure of the gripper is shown in Fig. 5a. The gripper 
is set in the fully open position and closes on a 2-cm diameter ball 
weighing 0.5 kg. The impact forces are quite high, several times the 
magnitude of the final static gripping force applied by the servomotors. 
Figure 5b shows what happens when the mass is reduced by positive 
acceleration feedback. The amount of positive feedback must be care- 
fully controlled or effective mass will become negative and result in 
oscillations. For the results of Fig. 5 the acceleration gain has been 
adjusted to reduce the effective mass to about 10 percent of the true 
mass. The impact forces have been reduced by acceleration feedback, 
but the amount of reduction is dependent upon properties of the 
rubber ball and other feedback parameters. Under good conditions, 
impact forces can be reduced by almost 50 percent. 

There are two significant differences in the responses of Figs. 5a 
and b. First, the initial peak impact force has been reduced from 50N 
to about 20N with acceleration feedback. Furthermore, the gripper 
stabilizes the ball position much more rapidly. This is due to the 
increased responsiveness of the gripper fingers to the ball dynamics 
as a result of the reduced effective mass. Increased motor drive 
bandwidth is also required in order to control the effective mass. Note 
the high frequency components present in the drive currents of Fig. 
5b. 

A series of test simulations were conducted for various rubber ball 
stiffnesses and damping coefficients to investigate the effectiveness of 
acceleration feedback. Again, simulations were performed using a step 
function position control input. The ball is positioned 1 cm off-center 
toward the right finger in order to test the ability of the gripper to 
center and stabilize the position of the ball. The position feedback 
gain was set to 500 amps per meter (which is equivalent to a spring 
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Fig. 5—(a) Damped gripper closure with position and velocity feedback. (b) Damped 
gripper closure with position, velocity, and acceleration feedback. 


rate of about 2750 N/m) for the first two test series, and critical 
damping was established by a velocity feedback of 14.95 A-s/m. 

The ball damping was set to 1.0 N-s/m for the first test series. The 
ball stiffness was varied from 10? N/m to 10° N/m. Acceleration 
feedback was set to reduce effective mass to about 5 percent of 
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mechanical mass at 0.11 A-s?/m. The ratio of peak impact force applied 
by the first contact between a gripper finger and the ball with and 
without acceleration feedback was taken to be the measure of the 
effectiveness of this feedback. The results of this test series are plotted 
in Fig. 6a. The plot shows that the ball stiffness has a substantial 
influence on the effectiveness of acceleration feedback. Peak impact 
force is only slightly reduced for very large or very small stiffnesses. 
The effect is substantial at moderate stiffness (10* N/m) but only over 
a narrow range. 

A second series of tests were conducted with the stiffness fixed at 
10* N/m. The ratio of impact forces was plotted again for ball damping 
from 0.01 to 10* N-s/m. Figure 6b shows that damping has little effect 
on the impact ratio until it exceeds about 10 N-s/m. After this point 
is reached, electrical damping predominates over the stiffness. At very 
large damping values the acceleration feedback becomes totally inef- 
fective. 

The position feedback was varied in a third set of test simulations. 
The ball stiffness was fixed at 10* N/m-and the damping was set at 1 
N-s/m. Velocity feedback was adjusted in each instance so that critical 
damping was maintained. Figure 7 shows a plot of the results for 
position feedback gains varying from 10 to 5000. Both very low and 
very high gains appear to adversely influence the effectiveness of 
acceleration feedback. Position feedback gains in the range of 50 to 
200 A/m have the greatest effect. 

Since position feedback produces an effect similar to the stiffness 
of the ball, it is not surprising that the plots have similar character- 
istics. Loss of acceleration feedback effectiveness at high stiffnesses 
can be attributed to relatively slow response time of the servomotors 
in comparison to the mechanical resonance of the ball and finger 
system. Since it takes approximately 100 us for the servomotors to 
reverse the gear force (because of backlash), the mechanism simply 
does not have enough time to counteract the effects of the real mass 
of the fingers. Poor acceleration feedback effectiveness at very low 
stiffnesses appears to be due to mismatch between the stiffnesses of 
the finger and the ball. If one of the stiffnesses overwhelms the other, 
particularly by several orders of magnitude, the acceleration feedback 
becomes useless. 

Another means for reducing impact forces is to bring the fingers 
into contact with the ball at low velocity. This could be accomplished 
by simply bringing the fingers slowly together by applying a ramp 
function to the position control inputs; however, this makes the gripper 
operation slow. Since it is only necessary for velocity to be low just 
before contact is made with the ball, another possibility is to increase 
the negative velocity feedback so that the finger system is overdamped. 
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Fig. 6—(a) Peak impact force ratio versus ball stiffness. (Cont.) 


This results in an exponential velocity response. However, computer 
simulation shows that this does not really reduce impact forces, since 
the fingers become very resistive to springy impacts. The series of Fig. 
6b illustrates this. Another alternative is to apply an exponential 
position control input to each finger. This takes advantage of the 
higher allowed velocity and slows the fingers down as they approach 
the contact points. Furthermore, these input functions are very easy 
to generate in hardware. Figure 8 shows an example of combined 
positive acceleration feedback and exponential control input functions 
with 50-ms time constants. The impact forces are now reduced to little 
more than the final static gripping force. Several time constants for 
the input exponentials were investigated, and the 50-ms exponential 
appeared to yield the best overall results (low impact force consistent 
with rapid closure) for this particular mechanical gripper. 

An analysis of sampling period was performed on a discrete time 
control architecture to determine if the gripper could be controlled 
adequately by a microprocessor. This was accomplished by modeling 
the motor drivers as zero-order sample and hold devices. A minimum 
sample rate of at least 5000 samples per second was indicated by 
numerous computed responses. This rate would be needed to control 


954 TECHNICAL JOURNAL, APRIL 1985 


0.9 


| 


0.8 


EFFECTIVE MASS 
REAL MASS 


0.7 


FORCE RATIO 


0.6 


0.5 
10-2 107! 1 10 102 103 


BALL DAMPING IN NEWTON SECONDS PER METER 
Fig. 6—(b) Peak impact force ratio versus ball damping. 


the gripper without loss in performance. This requirement could be 
relaxed somewhat if certain erroneous dynamic response characteris- 
tics could be tolerated; however, the characteristics observed in the 
simulation were not considered acceptable. This sample rate would 
only be possible for a microprocessor if a very simple control strategy 
were being implemented. Computing the motor currents necessary 
every 200 us for the more sophisticated feedback strategy being 
planned is more than a normal microprocessor can handle. (It is 
possible to perform the operations with a digital signal processor, 
however.) Consequently, the controller design is primarily analog, with 
digital-to-analog converters to allow a microprocessor to control the 
various feedback gains and apply position control inputs. 

A conference on robot manipulators’ triggered some thinking on the 
concept of differing common-mode and differential-mode compliance. 
Active compliance in our gripper is useful as long as the robot arm is 
not moving. When the robot arm does move, however, this previously 
useful compliance becomes a liability because the fingers of the gripper 
will swing from side to side under the acceleration influence of the 
arm motion. It is desirable to maintain the compliance relative to the 
grasped object and simultaneously keep the fingers and object from 
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Fig. 7—-Peak impact force ratio versus position feedback gain (critically damped). 


swinging. This can be accomplished with appropriate common-mode 
feedback as discussed in Section II. 

The effect of additional common-mode stiffness is illustrated in Fig. 
9. This should be compared to Fig. 5a. Both common-mode position 
and common-mode velocity feedback are necessary for suitable com- 
mon-mode performance. The initial conditions have been set to the 
same values as in previous tests. Inertia of the ball causes the gripper 
to accommodate somewhat to the position of the object and close 
slightly off-center. The object is then forced to the center by the 
common-mode stiffness. Close inspection of the position plot of Fig. 
5a shows some common-mode finger oscillation under the influence 
of the ball vibration. This common-mode component is almost entirely 
eliminated in Fig. 9. Further experiments simulating arm motion by 
applying common-mode variations to the position control input func- 
tions indicate that the increased common-mode stiffness significantly 
reduces the amount of arm motion influence on the fingers. 

We had previously assumed that velocity and acceleration signals 
were available from appropriate sensors such as tachometers and 
accelerometers. The actual gripper uses a position sensor that consists 
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Fig. 8—Damped closure with 50 milliseconds time constant. 


of a capacitive measurement device described by G. L. Miller.® The 
velocity and acceleration signals are obtained electronically from the 
position signal by electrical differentiation. It is desirable to limit the 
high-frequency gain of the differentiator circuits used to estimate 
these parameters, so electrical poles are purposely designed into the 
differentiator transfer functions. Extensive computer simulation has 
shown that the electrical time constants should not exceed 100 us. 
Figure 10 shows a time response for circuit time constants of 10 ms. 
This should be compared to Fig. 9. The obvious distortion of the 
mechanical response is further indicated in the log error plots of 
velocity and acceleration estimates for the right finger. The errors are 
typically no more than about 10 percent for circuit time constants of 
100 us, and the mechanical response is virtually identical to responses 
obtained with smaller time constants. 


V. TESTING THE GRIPPER 


Initial tests conducted on the gripper showed that high-frequency 
noise was present in the position signal. Although the level of this 
noise was low (about 50 dB s/n), it became important when used to 
determine the velocity and acceleration signals. At the acceleration 
output of the differentiators, the wideband noise level was only 30 dB 
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Fig. 9—Damped gripper closure with common-mode feedback. 


below the signal, a level high enough to cause gear chatter in the rack 
and pinion drive. Subsequent modifications of the position sensing 
circuitry reduced this noise level to about 90 dB below signal, which 
made the acceleration signal useful for feedback purposes. 

To verify the simulation results, measurements of the mechanical 
finger mass and damping were necessary. This was accomplished by 
performing underdamped step response tests on the gripper. Figure 11 
shows a typical underdamped step response. This response is very 
much like that of a simple mass, spring, and damper. The mass can 
be determined from the resonant frequency and the damping from the 
rate of decay in the response. Undamped resonant frequency is given 
by : 

pee (24) 
m 
where w is frequency in radians per second, k is the spring constant in 
Newtons per meter, and m is mass in kilograms. Damping (N-s/m) 
can be determined from 
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=== i, (25) 


where b is the damping in Newton-seconds per meter, m is the mass 
in kilograms, x; and x2 are successive peak amplitudes in the decaying 
response, and T' is the cycle period in seconds. 

There remains the problem of verifying the spring constant, which, 
in this case, is due to the position feedback through the servomotor. 
Although the motor is completely characterized by the manufacturer, 
it was deemed important to measure the spring constant actually 
obtained from the servomotor. This was accomplished by performing 
a step response analysis with two different masses attached to the 
motor through the rack and pinion, as follows: 


Let m, = finger mass 
my, = slide and motor mass. 
Then, 
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k 
ee 26 
oe mM, + Mo ( ) 
gives the undamped resonant frequency of the complete motor, slide 
bearing, and finger masses for a given equivalent spring constant R. 
The finger mass m, can be removed, leaving mz and a second resonant 
frequency w2. Thus an equation for k can be derived: 
my 


= (27) 


Sol 
S51 


Mass m, is measurable since it can be removed from the rest of the 
mechanism. Mass m, consists of a combination of the slide mechanism, 
the rack mass, and the gear and motor inertia. Once k is determined, 
then mz can be calculated from (26). 

Measurements in the laboratory yielded the following data: 


k = 3900 N/m 
m, = 33.12 g 
w, = 154.2 rad/s 
wo = 172.6 rad/s 
My, = 1381 g. 


A calculation of the gear inertia was made from the size of the gear 
and the density of steel (equivalent mass: 3.5 grams). Adding this to 
the published specification for the motor inertia and translating into 
the linear (rack) units yields a slide mass of 104 grams. Thus the total 
slide and finger mass is about 137 grams. The servomotor inertia with 
the gear attached is about 1.77 x 10° kg — m”. With these parameters 
determined, the simulation program was used to verify the resonant 
frequencies. 
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Simulation results compared closely with the actual performance of 
the gripper. Typical step response settling times of 30 ms were obtained 
for high-position feedback gains, as predicted by the computer simu- 
lations. Absolute positioning accuracy is dependent upon the position 
sensor accuracy and mechanical imperfections. The currently em- 
ployed sensor has a maximum error of about 0.5 mm. 

Several slight differences in response between the actual gripper 
and the simulation were noted. The computer simulation did not 
consider the effect of static friction on the dynamic response. In this 
case, the frictional force appears to be a constant independent of 
velocity. With viscous friction, the underdamped response decays 
exponentially. Static friction results in a linear decay. The actual 
gripper responses indicate that a slight amount of static friction is 
present in the mechanism. The amount of static friction present is 
negligible for normally damped responses (near critical damping) since 
the viscous friction is dominant under these conditions. However, 
static friction does become significant when the finger compliance is 
high. This results in a positioning error, since the force necessary to 
overcome the static friction is proportional to the spring constant 
applied to the finger. One well-known method of solving this problem 
is to introduce high-frequency dither into the motor drive. The ampli- 
tude of the dither must be just sufficient to exceed the force necessary 
to break the static friction free. If the dither frequency is sufficiently 
high, its effect will not be seen in the finger position due to the inertia 
of the mechanism. Experiments with dither in the laboratory indicate 
that static friction is not the only source of positioning error. 

Dither was found to be an effective means for obtaining higher 
positioning accuracy in the presence of static friction at low-position 
feedback gains (low spring constant); however, magnet detent action 
in the servomotor due to the magnetic reluctance of the armature 
poles also causes positioning error. The effect of this motor phenom- 
enon is the appearance of virtual potential energy wells at uniform 
intervals along the finger track. Thus, the positioning error can be 
positive or negative, depending on the commanded position and the 
relative position of the nearest potential energy well. This effect can 
be compensated by adjusting the command position slightly, but it is 
nearly impossible to position the finger at a potential energy peak 
since any slight change in position can cause a change in sign of the 
detent force, which requires dynamic compensation. The only practical 
solutions for this problem are to increase the position feedback (which 
provides the dynamic compensation) of the system or reduce the motor 
gear radius so that the detents are closer together. The difficulty in 
electronically compensating for the detent action arises from the fact 
that there is no way of distinguishing between internally and externally 
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generated finger displacements with the present sensors. The current 
version of the gripper exhibits positioning detents at about 4 mm 
intervals. 

The effect of acceleration feedback was tested by studying the 
underdamped step response of the gripper fingers with a known 
position feedback gain. The frequency of the response is an indication 
of the effective mass of the system. It was found that the resonant 
frequency could be increased by about 20 percent, indicating an 
effective mass change of 


= 0.64, (28) 
my, | 
or about a 36 percent decrease in mass. The simulation predicted a 
somewhat higher mass reduction (about 55 percent), but it had not 
considered the effects of position measurement noise, which limited 
stable feedback gain. 

Position measurement noise also becomes significant in the velocity 
feedback term when velocity gain is high (high damping). When 
position feedback gain is set to a high level, high-velocity gain is 
needed to obtain critical damping. Under these conditions the noise 
present in the position signal causes servomotor gear chatter. This 
gear chatter has no significant effect on the positioning accuracy or 
dynamic response of the gripper, but it does produce annoying sounds 
in the rack and pinion mechanism when the finger is not in contact 
with an object and is motionless. The high damping in the finger is 
not generally needed except when the finger is moving. Under com- 
puter control, it was found that the dynamic response of the finger is 
not affected if the high-velocity gain is present for only about 50 ms 
after the “move finger” command is initiated. Then the velocity gain 
can be reduced to a safe level for noise-free operation. 

To make control of the gripper easier for the user a set of computer 
routines were written in the C language. These procedures allow the 
user to set all feedback parameters directly or, at a higher level, to 
move the fingers without having to be concerned with determining the 
appropriate velocity feedback for critically damped response. Velocity 
gain for critical damping. can be determined from an equation of the 
form | 


b = kive — ko, (29) 


where b is velocity gain, c is position gain, and k, and k, are constants. 
Constant k, is due to the mechanical damping already present. Another 
slower procedure provides a method for opening and closing the fingers 
slowly by incrementally stepping the fingers to the command position. 
More sophisticated algorithms allow the user to specify a compliance 
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and gripping force to be applied when closing. The gripper is closed 
on the object and the applied force is increased to the specified amount. 


VI. SUMMARY 


A robot gripper and mechanical impedance control system has been 
constructed based on the results of computer simulations. The com- 
puter simulation is an accurate and useful tool for studying the 
dynamic behavior of the gripper. It allows one to investigate various 
performance properties of the mechanical system (i.e., conduct com- 
puter experiments) without taking the risk of damaging components 
of the system or incurring the expense of building experimental 
equipment. The simulation could be improved somewhat by including 
the effects of static friction, the servomotor detent action, and the 
control system generated noise. 

The controlled impedance gripper has an adjustable finger position- 
ing stiffness that can be set in steps of about 90 N/m up to 23,000 
N/m. Damping can be set accordingly by computer so that near 
critically damped response is obtained. The effective finger mass can 
be reduced by acceleration feedback to about 64 percent of its true 
value. This helps to reduce impact forces when the gripper fingers are 
closed on an object. 

A computer algorithm library was devised to provide a user interface 
to the gripper control system. With this library a user can establish 
operating parameters for the gripper and move the fingers with criti- 
cally damped responses without having to know the necessary velocity 
feedback parameters. Because of noise in the position measurement 
system, high-velocity feedback gains cause servomotor gear chatter. 
An algorithm has been developed to allow critically damped responses 
at high-velocity feedback gains by applying the velocity gain only 
during the transient period, after which the velocity feedback gain 
is reduced to a statically stable level. In this way the chatter was 
eliminated. 

Improvements can be made to the gripper and control system to 
reduce the effects of static friction and motor detent action. Dither 
(high-frequency perturbation) of sufficient amplitude can reduce the 
effect of static friction. Smaller gear diameters will allow more accurate 

positioning of the fingers since detent positions will be more closely 
- spaced in finger coordinates. 

Adaptive control strategies are likely to be useful in future studies 
of the problems of robotic gripping. This becomes particularly impor- 
tant for manipulating objects of varying masses and physical dimen- 
sions. Knowledge of the size of objects being grasped may not, in 
general, be known in advance. This means that the gripper will have 
to adjust position feedback and/or control inputs to maintain a desired 
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gripping pressure on an unknown object. Further adaptations may 
improve the ability of the gripper to maintain the object in a stable 
position as the robot arm is moved, even if the mass of the object is 
unknown. In the future these operations will probably be performed 
by a microprocessor because of the relatively sophisticated techniques 
involved, but the analog controller designed here should provide the 
microprocessor with the control interface necessary to generate the 
appropriate gripper response. 
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APPENDIX A 
Runge-Kutta Method for the Solution of Differential Equations 


Runge-Kutta methods for solving general differential equations are 
well known and probably the most widely used algorithms for this 
purpose. The particular algorithm employed in this study is a fourth- 
order formula with coefficients attributed to Kutta (see Ref. 9). The 
algorithm solves the basic first-order equation: 


x(t) = f(x, t) (30) 


for x(t), where t is the independent variable. Since any high-order 
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differential equation can be written as a set of first-order differential 
equations, an equation of arbitrary order may be solved by repetitive 
application of the first-order algorithm. 

Equation (1) is integrated over a small interval h of the independent 
variable, which is termed the integration interval or step size. Each 
interval starts at an initial condition (e.g., x(0)) and concludes at an 
estimated value for x(t) at the end of the interval, [e.g., x(h)], which 
provides the initial conditions for the next integration interval. This 
process is repeated, yielding x(kh) for k = 1, 2, --- ,n. For systems of 
first-order differential equations, the integration is applied to each 
equation in parallel for each step. 

The fourth-order integration formula used is 


(t+ h) = £(t) + h(Ry + Qkp + hs + ky)/6 (31) 
ky = f(x, t) (32) 
ko = f(x + hki/2, t + h/2) (33) 
hs = f(x + hkp/2, t + h/2) (34) 
ky = f(x + hks, t +h). (35) 


Four evaluations of f(x, t) are required for each integration step. This 
formula reduces to Simpson’s rule if f(x, ¢) is a function of t only. 

The integration truncation error committed at each integration 
interval is dependent upon the size of the interval h and for an mth- 
order Runge-Kutta method the error is, in general, O(h™*'). The 
actual error committed is further dependent in a complicated manner 
upon the dynamics of the equation. More specifically, a smaller inter- 
val will be required to maintain the same error bounds during rapidly 
varying periods of the solution function x(t) than will be required 
when the solution function is slowly varying. Since f(x, t) must be 
evaluated four times for each integration interval and f(x, t) may be 
quite complicated and computationally expensive, it is desirable to 
keep h as large as possible without incurring excessive error accumu- 
lation. For this reason the integration interval is adaptively increased 
and decreased as the integration proceeds. 

To determine the proper integration interval during the integration 
process, an estimate of the truncation error being committed at each 
step is necessary. Several methods for estimating truncation error and 
adjusting step size have been reported in the literature (e.g., Collatz’®); 
however, these methods are often nearly as time-consuming as the 
integration itself or suffer from certain serious failure modes. The 
Collatz method, for example, is only a qualitative measure, which is 
quite efficient, but suffers from a divide-by-zero problem when the 
slope of x(t) becomes zero. A new method employed in this implemen- 
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tation of the algorithm was developed by the author to avoid excessive 
computation and eliminate the possibility of failure due to special 
cases of the solution. This method is also a qualitative measure and 
does not give a very accurate estimate of the true truncation error, but 
empirical study indicates it does yield a figure of merit that is within 
about an order of magnitude of the true error. By specifying figure-of- 
merit bounds that are an order of magnitude smaller than the maxi- 
mum error required, we can easily maintain a sufficiently small inte- 
gration interval without excessive computation or other complications. 

This error estimator is based on a measure of the difference in 
truncation error expected between a second-order Runge-Kutta 
method and the fourth-order Runge-Kutta. Equation (35) is a second- 
order predictor for the solution at the end of the integration interval, 
x(t + h). Recall that the truncation error committed at each step by 
an mth-order integration is O(h”*'). By comparing k, to the solution 
obtained from the fourth-order integration and scaling the result by 
h?, an error estimate proportional to the fourth-order truncation error 
is obtained. That is, 


error = a(h’)|k, — x(t + h) |, (36) 


where x(t + h) is the fourth-order solution for f(x, t), h is the 
integration interval, and a is a constant (usually between 0.1 and 10). 
If the error estimate exceeds a predetermined bound, then the integra- 
tion interval is reduced to one-half its current value. Likewise, if the 
error estimate is more than (say) a thousand times less than the error 
bound, then the interval is increased to twice its previous value. 


APPENDIX B 
State Equation for the Two-Fingered Gripper 


The following listing is written in the C programming language. It 
describes the complete state equation for the two-fingered gripper with 
common-mode feedback and the electrical time constants. 


int gripper (t, y, dy) 
doublet, y[], dy[]; 
{ a 
double u1(), u2(); 
double fabs(), 10910(); 
double kg, pg, f, f1; | 
_/* state variables for the robot gripper: 


0. left finger position 


1. left finger velocity 
2. left motor angle 
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3. left motor angular velocity 
4H. right finger position 
5. right finger velocity 
6. right motor angle 
7. right motor angular velocity 
8. rubber ball position 
9. rubber ball velocity 
10. left vel. RC network 
11. left acc. RC network 
12. right vel. RC network 
13. right acc. RC network 
14. common-mode vel. RC network 
15. common-mode acc. RC network 
16. left motor driver RC 
17. right motor driver RC 
* / 
kg = 4.16e6; /* spring constant and damping coefficient 
for gears */ 
pg = 322.32; 
/* compute state derivative function for left hand finger 
system */ 
dy [0] =y[1]; 
if(fabs((f=rg* y[2]—y[0]))>H){ /*gear force */ 
Fix=(f£>0. ? (f-H) : (f£+H)) *kg+ (rg * y[3] —y[1]) * 
Pg; 
if(f>0.){ 
fi=(fi>0. 2? £1: 0.); 
lelse{ 
fi=(f£1<0.?£1:0.); 
lelse{ 
f1=0.; 


, 


if((if=y[0]—y[8]+R)>0.){ /* rubber ball force */ 
£2=(f *k3+(y[1] —y[9]) * P3); 
f£2=(f2>0.?fs:0.): 

lelse{ 
£2=0.; 

dy[1]=(-P * y[1] +£1—£2)/M; 

/* compute left motor current and voltage */ 
Iai=y[16]; 
if(fabs(Ia1) >Imax) { 
Iai=(Ia1>0. ? Imax : —Imax) ; 
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V1=Ra * Ia1+La * (Iai —Iaisv)/dt+Ka * y[3]; 
if (fabs(V1) >Vmax) { 
Vi=(V1>0. ? Vmax : —Vmax); 
Iai=(v1+La * Iaisv/dt —Ka * y[3]/(Ra+La/dt); 
dy[2] =y[3]; 
dy[3] = (Km* Ia1—p*y[3]—f£1*rg)/m; /*torque */ 
/* compute state derivative for right hand finger system */ 
dy[4]=y[5] 
if (fabs ((f=rg* y[6]—y[4]))>H){ /* gear force */ 
fl=(f>0.? (f—-H): (E+H)) *kgt+t(rg*y[7]—y[5]) * 
Pg; 
if (£>0.){ 
fi=(f1>0.?£f1:0.); 
lelse{ 
fi=(f£1<0. ? £1: 0.); 
lelse{ 
fi=0.; 
if ((f=y[8]+R-—y[4])>0.){ /*rubber ball force */ 
f£3=f£ *k3+ (y[9]—y[5]) * p3; 
f3=(f3>0.?7?f3:0.); 
else f{ 
£3=0.; 
| 
dy[5J=(-P * y{[5]+£1+£3)/mM; 
/* compute right motor current and voltage */ 
Ia2=y[17]; 
if (fabs(Ia2) >Imax) { 
Ia2=(Ia2>0. ? Imax : —Imax); 
V2 =Ra * Ia2+La * (Ia2 —Ia2sv)/dt+Ka * y[7]; 
if(fabs(V2) > Vmax) { 
V2=(V2>0. ? Vmax : —Vmax); 
Ia2 = (V2+La * Ia2sv/dt —Ka * y[7]/(Ra+La/dt); 
dy[6] =y[7]; 
dy[7] = (Km* Ia2—p* y[7] —£1* rg)/m; 
/* compute state derivative for rubber ball */ 
dy[8]=y[9]; 
dy[9] = (£2 —£3)/m3; 
/* state equations for electronic circuitry */ 
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dy[10] =(y[0] —y[10])/Rc; 

dy[11] = (dy[10] —y[11]/Rc; 

dy[12] =(y[4] —y[12])/Rc; 

dy[13] = (dy[12] —y[13])/Rc; 

dy{14] =((y[0]+y[4]) —y[14])/Rc; 

dy[15] = (dy[14] —y[15])/Rc; 

dy[16]=(-a *dy[11]-—b*dy[10]+(c+£) * (u1(t)—y[0]) 
—d*dy[15]—e*dy[14] +f * (u2(t) —y[4]) —y[16])/ 
1.2e—4; 

dy[17] = (-a * dy[13] —b * dy[12] + (c + £) * (u2(t) -—y 
[4]) —-d * dy[15] —e * dy{[14] + £ * (ui(t) — y[0]) —- 
y[17])/1.2e—-4; 

return (18); 
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Determination of Fiber Proof-Test Stress for 
Undersea Lightguide Cable 


By T. C. CHU and H. C. CHANDAN* 
(Manuscript received July 24, 1984) 


A method is developed to determine the fiber strength requirements for the 
undersea lightguide cable. This method utilizes the theories involved in fiber 
proof testing, undersea cable dynamics, and nonlinear tensile behavior of 
cable. The method is illustrated with a sample cable currently being developed 
for the transatlantic system. This method provides great flexibility in design- 
ing undersea lightguide cable. It can be applied to justify a new cable design 
at an early stage or to select proper fiber proof-test levels for an existing cable 
under different operating conditions, i.e., deeper ocean, worse sea state, or 
faster recovery. 


I. INTRODUCTION 


Because of the difficulty and high cost of repair, undersea cable 
requires high reliability. Operations such as laying, recovery, and 
holding impose high tensions on the cable.’ The high tension causes 
high strain and permanent deformation in the cable, and it may reduce 
or impair system reliability. Much experience has been gained over 
the past decades in designing conventional undersea cable with me- 
tallic conductors to withstand the high tension. In general, cable 
strains well in excess of 1 percent can be tolerated for the coaxial 
cables with metallic conductors due to the ductile nature of metals. 

For the new undersea cable that has lightguides as the signal- 
carrying medium, the problems of high tension and strain cause greater 


* Authors are employees of AT&T Bell Laboratories. 
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concerns because of the static and dynamic fatigue of optical fibers.”* 
To ensure that the fibers survive the high-tension cable operations 
(laying, recovering, and holding) and last for 25 years, the entire length 
of each fiber is proof tested” to a given stress before cabling. The 
proof-test level is chosen by taking into account the load-time history 
of the fibers during the entire service period. An appropriate proof- 
test level must be chosen since an underestimated proof-test stress 
reduces the system reliability, and an overestimated one increases 
system cost due to reduced fiber yield or increased fiber loss caused 
by more frequent repair splices during proof testing. 

In this paper a method is developed to determine the fiber proof- 
test stress for undersea lightguide cable. In Section II the theory of 
proof testing is discussed, with special attention to the material 
constants used in the calculations. In Section III the cable tensions as 
function of time are analyzed. In Section IV the nonlinear tensile 
property of the cable is discussed and is used to obtain the cable and 
fiber strain. The fiber strain is then used to determine the fiber proof- 
test stress. In Section V the method is applied to a sample cable that 
is currently under development for transatlantic systems. 


Il. PROOF-TEST LEVEL 


It is well established that the strength of glass optical fibers is 
determined by micrometer- or submicrometer-size surface or subsur- 
face flaws. These flaws grow with time under the combined action of 
tensile stress, temperature, and ambient moisture. These fracture- 
initiating flaws are randomly distributed along the length of the fiber. 
Therefore, the entire length of the fiber is proof tested to screen out 
flaws larger than those that correspond to the proof-test level. For 
example, the largest initial flaw present in a given length of 200 kpsi, 
proof-tested, fused-silica fiber is 0.212 um.* 

Subsequently, the proof-tested fiber is subjected to varying stresses 
for different time periods during cable manufacture, laying, recovery, 
and repair. During these operations, the preexisting flaws in the fiber 
grow, and this fiaw growth is cumulative. At present, no commercial 
scheme is available that will prevent crack growth in fibers. Therefore, 
we choose the proof-test level that restricts crack growth during 
expected worst-service conditions to yield a desired design life for the 
fibers. In other words, by proof testing we ensure that the initial crack 
size is sufficiently small so that the full range of stresses to which the 
fiber is subjected in service will not cause the fiber to rupture. 

Assuming a power law for crack growth, the largest crack, a,, in a 
proof-tested fiber will grow at a rate (da)/(dt) given by” 

da 


Ga _ AKN 
dt AKj ’ (1) 
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where K; is the stress-intensity factor (the subscript J refers to the 
crack-opening mode) at the crack tip and is given by 


Ky = Yoa"?, (2) 


where Y is a geometrical factor related to the shape of the crack and 
the specimen (~1.241 for semicircular part-through crack, normal to 
the axis of a cylindrical fiber); o is the applied stress. The characteristic 
crack-growth parameters A and N are coupled and are empirically 
determined from dynamic and static fatigue tests.° From a collection 
of experimental data covering several independent investigations, 
we have observed that the relationships between log A and N are 
different for dynamic and static fatigue situations. For short-term 
loads Ag, increasing with time, our studies have indicated that log A, = 
3.22N — 8.516. For a static load A,, log A, = 3.289N — 10.05. The 
choice of N depends on the Relative Humidity (RH) in the ambient 
environment: N = 14 to 16 for wet environment (97-percent RH) and 
N = 17 to 20 for typical laboratory environment (40- to 60-percent 
RH). 

From eqs. (1) and (2), an applied stress of o; for a time ¢, will cause 
the crack to grow from the proof-test level of a, to a, and can be 
calculated by 


N N ty 
a? ay? ==T> | Ay [a(t)]Nat, (3) 
0 


where a, is the flaw size after the fiber has been under an applied 
stress o;(t) for a time ¢t,. The fused silica optical fibers behave in a 
linear elastic manner up to the breaking point; therefore, the stress 
o,(t) in eq. (8) can be expressed in terms of strain, ¢,(t) = [o1(t)]/E, 
where E is the Young’s elastic modulus of the fiber (E = 71.9 GPa or 
10.4 X 10° psi). Furthermore, when the fibers are subjected to a 
sequence of stress loadings, the flaw growths are additive.* From these 
considerations, eq. (3) can be rewritten as 


Q ti 
ap = af — ¥ amexy® | [e:(t)]" dt, (4) 


i=1 


where M = (2 — N)/2 and aq; is the final flaw size after the sequence 
of stress loadings. If the flaw growth is considerable, say a;/a, = 10, 
and we assume that JN is in the range of 14 to 20 (M = —7), then 

a; 

— 72107. 
ap 
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Thus, the term a/ can be neglected and 


a E>) 


M ~ 
M = 


a 


AMEN YY |" [a(n (5) 


i=1 


The required proof-test level is determined from the fracture mechan- 
ics relation, : 


Kic 
Op — Ya}? ’ (6) 


where Kic is the fracture toughness of the material. For fused silica 
fibers, 


Kic = 0.789 (unit in MPa Vm), 
and as previously mentioned, Y = 1.241. Hence, 


0.6357 
op — que ? (7) 
p 





where a, is determined from eq. (5) in meter and o, is in MPa. Thus, 
the fibers proof tested at o, will endure the stress loadings during 
cable laying, recovery, repair, and long-term service without fracture. 


Ill. CABLE TENSIONS DURING THE OPERATIONS 


To evaluate the fiber proof-test stress, the fiber strain ¢ in eq. (5) 
has to be determined. The amount of fiber strain is related to the cable 
tension, cable tensile property—i.e., tension versus strain relation— 
and the coupling between the fibers to the cable structure. In a cable 
with its fibers fully coupled to the cable structure, the fiber strain is 
equal to the cable strain, and this is the worst case. By knowing the 
cable tensions as functions of time during each operation, the cable 
strain as a function of time can be obtained from the cable tensile 
property. The cable tension as a function of time can be derived from 
the maximum tension at the ship, as shown below. 


3.1 Laying 

In a normal cable-laying operation, the cable forms a straight-line 
configuration, as shown in Fig. la. The cable tension at the ocean 
bottom is theoretically equal to zero. The cable tension at the ship is 
equal to wh,’ where w is the cable weight per unit length in seawater 
and h is the ocean depth. By following a cable element during the 
laying operation, the time-dependent tension experienced by the ele- 
ment can be found as 
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he E 


0<t=<=——,, 
h vsin(as) 
vsin(a,) 


where v, is the cable-laying speed, t is the elapsed time after the cable 
element leaves the ship and enters the water, and a, is the angle 
between the cable and sea level and is 

H 


Aas = ~s, 


U! 


(8) 


where H is the hydrodynamic constant of the cable. For most undersea 
cable used in deep water, H ranges between 30 to 50 degree-knots. The 
cable-laying speed is usually about 7 to 8 knots, thus, the angle a, is 
less than 10 degrees. For such small angles, the contributions of wave 
motion and repeater weight to the tension are negligible. Equation (8) 
is sketched in Fig. 2a. 








(b) 


(c) 
Fig. 1—Cable operations. (a) Laying. (b) Testing or splicing. (c) Recovery. 
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T(t) 








wh 
Be re prea eee 
eo h 


vp sin (a) 


salt 2. 
Vg sin(a,) 
T(t) 


(b) 


a 





5. 
U, 


Fig. 2—(a) The time-dependent tensions experienced by the cable element during 
the laying operation. (b) Upper bound for cable tension during testing. (c) Catenary 
formed by the cable during recovery. 


3.2 Recovery 


During recovery, the cable forms a catenary, as shown in Fig. 1b, 
and the length of the catenary S can be evaluated.’ For a given recovery 
speed v,, the recovery time of a cable element is given by 


S 
t=—. 


Ur 


Since the cable forms a catenary, the cable tension at the ocean 
bottom, TJ», is not zero. The cable tension at the ship, not including 
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the effects of wave motion and repeater weight, is Mwh, where Misa 
magnification factor that indicates the difference in tension between 
laying and recovery. The factor M has the following form:' 


1 


M = er 
af cosa+cosa, \= 
1 — tan? — SR ra AY 
2 \1 — cos a cos a, 
sin’ a 
Of. Se ee 
2 —sin?a’ 


where a, is the angle between the cable and sea level and a is the ratio 
between the cable hydrodynamic constant H and the recovery speed 
v,. Assuming that the tension experienced by a cable element during 
recovery varies linearly with time, the cable tension as a function of 
time—including the wave motion and repeater weight—turns out to 


be 


T(t) = Ty + ne as ey le 
U, 


o<t<, (9) 


where T;, is the weight of a cable-carried point mass (such as a repeater) 
in seawater, and T,, is the additional tension due to the wave motion. 
T,, is expressed by 


—— 2rt 
Ty = VEAp,vyCOS = COS Qs, 


where EA is the cable tensile stiffness, p, is the cable mass per unit 
length, v, is the maximum ship vertical velocity, 7 is the period of 
wave motion, and a, is the recovery angle. In normal recovery, the 
angle a, is kept close to 90 degrees to reduce the tension, thus cos 
a; = 1. Equation (9) is expressed in Fig. 2b. 


3.3 Holding 


After its recovery, the cable is held for a period of time for testing 
and repair. At present, the required splicing time for the lightguide 
cable is estimated at 16 hours. Assuming (pessimistically) that three 
attempts are made before a successful repair is executed, then the 
total holding time is 48 hours. During holding period, the cable also 
forms a catenary that 1s essentially stationary except for the effect of 
wave motion, as shown in Fig. Ic. The tension experienced by the 
cable element at the ship can be shown as 
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T(t)=T + Ty 0<t< bt, (10) 
where fy is the holding time, T,, is the additional tension due to the 
wave motion, 7’; is the tension due to the catenary and can be evaluated 
from the theory of a catenary as’ 


T,. = wC cos h(z2), 


where w is the cable weight per unit length in water, and C and z are 
constants to be determined from the ocean depth h and the cable angle 
a; at holding. Because the cable sections suspended near the ocean 
bottom experience lower tension than the cable elements near to the 
ship while holding, the survival of cable elements at the ship guaran- 
tees the survival of the remaining cable. In short, T(t) = T, + T, 
expresses the upper bound for cable tension while loading. This 
equation is represented in Fig. 2c. 


IV. CABLE TENSILE PROPERTIES 


By knowing the cable tensions during the operations, the cable 
strain can be found if the tensile property of the cable is known. A 
method has been developed to accurately predict the tensile property 
of undersea cable in both elastic and plastic regions.® The method 
relates the cable tension T to cable strain e, by a fifth-degree poly- 
nomial as 


T = Co + Cre, + Coe? + Cge2 + Cred + Cre2, 


where Cy through Cg, are constants depending on the mechanical 
properties of the constituent cable components. An example of such a 
polynominal is shown in Fig. 3. Since our present need is to find the 
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Fig. 3—Tensile properties of the TAT-8 transatlantic cable. 
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Fig. 4—Cable structure and dimensions for the TAT-8 transatlantic cable. 


cable strain for a given tension, an inverse function of the above 
relation is required. This can be done by the numerical regression 
method. The resulting relation is 


re Oe ea Ord Ae om Ory Gee errr oa Get Aus (11) 


where k and Cy through C;, are determined with desired accuracy from 
the numerical regression. 

Since the cable tension during its operations is a function of time 
T(t), eq. (11) gives the cable strain as a function of time, i.e., «(t). If 
the fibers are fully coupled to the cable structure, the fiber strain is 
equal to the cable strain «(t) = «,(¢), and this is the strain required in 
eq. (5) to evaluate the fiber proof-test level. 


V. APPLICATIONS 


This method of testing fiber strength is applied to the cable that is 
currently being developed for the TAT-8 transatlantic system. The 
cable construction and dimensions are shown in Fig. 4. The cable and 
fiber properties required in the proof-test stress evaluation are listed 
in Table I. The cable tension during each operation is evaluated under 
realistic operating conditions. 


Table |—Cable and fiber properties 





Cable weight per unit length in seawater w 0.326 lb/ft 
Cable hydrodynamic constant H 43.5 degree-knots 
Cable mass per unit length p, 0.0174 —— = 
Cable tensile stiffness HA 2.9 x 10° Ib 
Fiber “N” number 17 


Fiber tensile modulus I 10.4 x 10° psi 


UNDERSEA LIGHTGUIDE 979 


5.1 Laying 


Following are the conditions under which cable was laid: 

1. Laying speed v; = 7 knots 

2. Ocean depth h = 18,000 ft (transatlantic installation) 

3. Cable angle a, = H/v; = 6.22 degrees 

4, Cable touchdown time = h/[vsin(a,)] = 3.93 hours. 
Considering the above conditions and the cable properties listed in 
Table I, eq. (8) becomes 


5875 
T(t) = 9875 — 303 x 60 x 60 
0 < t < 3.90 X 60 X 60 seconds. (12) 


5.2 Recovery 


Following are the conditions under which the cable was recovered: 
1. Recovery speed v, = 1 knot 
2. Recovery angle a, = 85 degrees 
3. Vertical ship velocity due to wave motion v, = 17 ft/s (sea 
state 6) 
4. Ocean depth h = 18,000 ft 
5. wh = 5875 lb 
6. Additional tension due to wave motion T,, = 3829 lb. (To 
simplify the integration, the tension due to the wave motion is assumed 
to be equal to V FAp.v,, the amplitude of the oscillating tension. This 
simplification results in a more conservative fiber proof-test stress.) 
7. Mwh = 138,488 lb 
8. Cable bottom tension Ty = Mwh — wh = 7613 lb 
9. Suspended cable length! = 1.65 x 18,000 = 29,700 ft 
10. Recovery time t = S/v, = 4.89 hours 
11. Repeater weight in seawater T= 1000 lb (presently, the repeater 
weight has been reduced to 300 lb in water). 
With the above conditions and cable properties, eq. (9) becomes 


T(t) = (7613 + 3829 + 1000) + 0.33373t 
0 < t <= 4.89 x 60 X 60 seconds. (13) 


5.3 Holding 


Following are the conditions during holding: 

1. Ocean depth h = 18,000 ft 

2. Cable angle at holding a, = 75 degrees 

3. Additional tension due to wave motion T;, = 3829 Ib 

4, Cable weight per unit length in water = 0.326 lb/ft. 
By knowing the ocean depth h and cable angle a,, the constant C and 
z can be evaluated from the catenary theory as 
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C = 6299.5 ft 
z = 2.0563 
T. = WC cos h(z) = 8133 lb. 
Thus, eq. (10) becomes 
T(t) = 8133 + 3829 
0 <t < 48 X 60 X 60 seconds. (14) 


5.4 Cable tensile properties 


The theoretical tensile property of the submarine cable evaluated 
by the method previously discussed is shown in Fig. 4.2 The inverse 
polynominal is found to have the following form: 


e= Co + Cif + C3T* + C,T* (15) 
for 0 < T <= 21,300 lb, where 
Co = —2.9758111668 x 10° 


C, = +3.96935593146 x 107’ 
C, = —8.1189437448 x 107” 
3 = +6.40675624028 x 10° 


and the standard error x = 6.9 X 10~°. It is noted that eq. (15) applies 
for cable tension below 21,300 lb. It is adequate because the maximum 
tension expected by the cable during its operations is about 17,500 lb. 


5.5 Proof-test stress evaluation 


The cable tension 7'(t) in each of eqs. (12), (13), and (14) is 
substituted into eq. (15) to obtain e,(t), e2(t), and e3(t). Substitution 
of €,(t), e2(t), and e3(t) into eq. (6) for «;(t) gives, upon numerical 
integration, the initial flaw size a,: 

az’ = 5.5685763 X 10° + 5.139996 x 10* + 5.0337450 x 10” 

(laying) (recovery) (holding) 
or 


A, = 3.2 X 10°-’meter = 0.32 um. 
The corresponding proof-test stress is obtained from eq. (7): 
Op = 1.124 GPa = 163 kpsi. 


Thus, a proof-test stress greater than 163 kpsi should guarantee the 
survival of the fibers in the cable. Generally, however, a somewhat 
higher proof-test level (200 kpsi) is specified for transatlantic cable to 
ensure that adequate margin is provided to offset any subtle effects 
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brought about by the cable manufacturing and handling processes that 
are not covered by the model. 


VI. CONCLUSIONS 


A method is developed to evaluate the required proof-test stress for 
the fibers used in the undersea cable. The method provides great 
flexibility in designing the cable. It can be applied to justify a cable 
design and verify the adequacy of a fiber proof-test level at an early 
stage, or to select proper fiber proof-test levels for an existing cable 
under different operating conditions, i.e., deeper ocean, worse sea 
state, or faster recovery. 
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Beamwidth and Useable Bandwidth of Delay- 
Steered Microphone Arrays 


By J. L. FLANAGAN* 
(Manuscript received October 31, 1984) 


Automatic delay steering of microphone arrays improves sound pickup in 
large rooms. But steering to wave-arrival directions away from broadside 
degrades the acuity of the beam and diminishes the useable bandwidth of the 
array. This paper derives quantitative relations for the variation in beamwidth 
and bandwidth for one-dimensional, uniform, unweighted arrays composed of 
(2N + 1) receivers spaced by distance d. The results show how the upper and 
lower useful frequencies and the beam acuity are conditioned not only by 
receiver spacing and frequency, but also by wave-arrival direction and steering 
direction. The relations developed permit detailed design of steerable arrays 
for specified frequency range and spatial coverage. 


I. MICROPROCESSOR CONTROL OF MICROPHONE ARRAYS 


Speech transduction or sound pickup in large rooms—such as au- 
ditoria or classrooms—traditionally is plagued by the distortions of 
room reverberation and noise interference from unwanted sources. 
The problem is minimized by obtaining, at the transducer, the greatest 
intensity possible for the signal travelling the direct path from source 
to receiver, and the least intensity possible for multipath room reflec- 
tions of that signal and any additive interfering noise. Arrays of 
microphones, especially two-dimensional arrays designed for high 
directivity, are useful for achieving this high ratio of direct to distorted 
sound. 


* AT&T Bell Laboratories. 
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But in an auditorium or conference room the desired sound source 
may shift position. The directional array must be steered, or pointed, 
to the shifting source to gain the full benefit. Rapid, automatic, 
electronic steering of the array is therefore attractive. Currently 
emerging microprocessors have the computational capability to per- 
form beam forming and beam steering dynamically, and can be pro- 
grammed to seek and track talkers in a room.’ 

Beam steering is implemented conveniently by digital delay adjust- 
ment of the signals received at each microphone of the array. But 
forming a beam and steering it to angles away from the normal to the 
array (i.e., away from “broadside”) have two effects that must be taken 
into account when designing the system. One is that the useful 
bandwidth diminishes as the steering direction is made acute to the 
axis of the array (or to the plane of the array for a two-dimensional 
system). Another is that the beam pattern becomes asymmetric and 
less sharp. This report considers the variation in beamwidth and 
useable bandwidth under conditions of delay steering. 


Hl. TIME, FREQUENCY, AND SPATIAL RESPONSES FOR THE UNIFORM 
DELAY-STEERED ARRAY 


The most convenient expository vehicle is the one-dimensional 
uniform array, unweighted in amplitude. Related analyses apply for 
two-dimensional and three-dimensional arrays. 

Figure 1 shows a line of (2N + 1) receivers, spaced uniformly by 
distance d. All receivers have the same sensitivity, taken to be real 
and equal to unity. For beam forming and steering, the individual 
receiver outputs are passed through controllable delays and are 
summed to produce the array output. 

For a plane, impulsive sound wave of unit amplitude, arriving from 
the direction of polar angle ¢, the time-domain and frequency-domain 
responses of the array are, respectively, 


N 
h(t)= ) 6(t+nT), 


n=—N 
and 
N . 
H(jw)= Y en, (1) 
n=—-N 
where 
t is time, 


w = 2rf is the radian frequency, 

6(-) is the delta function, and 

T = (7 — 7’) is the time difference between the interelement wave 
transit delay and the interelement steering delay, in which 
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T= - (cos¢) and 7’= : (cos 9’), 


where ¢’ is the steering angle, and c is the sound velocity (taken here 
as 3.4 X 10‘ cm/s). For notational convenience, a constant delay term 
that ensures formal causality has been omitted from (1). Further 
deserving of emphasis, the impulse response h(t) and its Fourier 
transform H(jw) are both implicit functions of the wave-arrival direc- 
tion ¢ and of the steering angle ¢’. For specific frequencies w, the 
spatial directivity is conveniently revealed in H(jw). An elaboration 
of the expression for H(jw) is therefore useful for our purposes. 
The finite geometric series of (1) can be reformulated as 
N 


H(jo) = ein 


n=—N 


(ef NTT 


= e JeNT) 


(e/°? — 1) 


sin[(2N + 1)(wT/2)] 


7 sin(wT/2) (2) 


WAVE-ARRIVAL 
DIRECTION 





o————-—0 


Fig. 1—Geometry for a one-dimensional microphone array of (2N + 1) elements. The 
microphone spacing is d, and the polar coordinate for wave-arrival direction is ¢. The 
array output is a summation of the outputs of individual microphones after delay 
adjustment. 
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2N + 1 


H(ja) 


Ty [AAAS 


wT/2 


Fig.2—Output of the uniform array characterized by sin[(2N + 1)(w7/2)]/ 
sin(wT'/2). 


Given values of the spatial coordinate ¢ and the steering parameter 
o’, one notices that H(jw) is periodic in w with period equal to 27/T. 
This relation is illustrated in Fig. 2. Low amplitude values of H( jw) 
represent conditions where the array is discriminating against sources 
in direction ¢ and for frequency w. At very low frequencies, by virtue 
of the characteristics of (2), the array cannot discriminate as a function 
of wave-arrival direction. At high frequencies, near 1/T and multiples 
thereof, a similar lack of spatial discrimination must occur. This 
“spatial aliasing” of the discriminatory capabilities fixes an upper limit 
to the useful frequency bandwidth of the array. This upper limit is 


Tissee = 1/| Tags 


where 
d / 
| T'| max = ‘ |cos @ — cos @’ | max. (3) 


(Unlike time-domain sampling, where signal and impulse train are 
multiplied and the signal spectrum translated to sidebands about the 
sampling frequency and its multiples, this spatial aliasing acts only as 
a filter for the original frequency components of the source signal.) 
The first zero of the response H(jw) can be taken as a measure of 
the lowest frequency below which the array cannot provide spatial 
discrimination. This first zero, as seen from (2), corresponds to 


(2N + 1)(wT/2) = x, 
or 
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_ __ hopper 
hower a (2N i 1) ‘ (4) 


As apparent from (3), given the directions ¢, ¢’, the upper useful 
frequency of the array is conditioned by the receiver spacing d. 
Similarly from (4), the lowest useful frequency for spatial discrimina- 
tion is dictated by the array size N. Alternatively, (4) gives the required 


N as 
1 fupper ) 
N=|= (2 —-1)) , 
E [Be nhi 


where nhi equals the next highest integer. 

Further, using the relevant zeros of (2) (i.e., the first, and the first 
less than fupper), AN approximation of the useful discrimination band- 
width of the array is 


BW = ( Fuppee — 2 fiower) + (5) 


Additionally, one notes that the value of | T'| is conditioned by the 
values ¢, ¢’ and has the limits 


0<|T] < 2d/c, 
corresponding, respectively, to 
@ = ¢’ (i.e., beam steered to the wave-arrival direction), 


and to 


@=0° for @’ = 180°, or 
@ = 180° for ¢’ =0° (6) 


(i.e., the beam steered to one axial direction with the wave arrival 
from the opposite axial direction). As a consequence of the trigono- 
metric relations, the useful upper and lower frequencies do not vary 
linearly with ¢’. This variation is plotted in Fig. 3 for a value of N = 
10. 

Finally, if w and @’ are prescribed as parameters, then (2) gives the 
amplitude response of the array as a function of the wave-arrival 
direction ¢. The response can then be considered an explicit function 
of the spatial coordinate, or H(¢). Wave arrival from the steering 
direction, ¢ = ¢’, gives an infinite period to (2). Also for this condition, 
H(¢) = | H(@) |max = (2N + 1), and fupper — ©, and the array passes 
all frequencies ideally. Wave arrival from a direction away from the 
steering direction diminishes the period of (2), as well as the output 
amplitude of H(@) for frequencies between the upper and lower values. 
This diminution in output reflects the beam discrimination, or beam- 
width, of the array. The half-power output (ie, |H(¢)|-3sa = 
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FREQUENCY IN HERTZ + c/d 


USEFUL 
BANDWIDTH 


7 flower (NV = 10) 
7 





90 60 30 0 
BROADSIDE ; END-ON 
STEERING ANGLE @ IN DEGREES 


Fig. 3—Variation of the useable array bandwidth as a function of the steering 
direction @’. 


0.707|H(¢ = ¢’)| is frequently used to define the beamwidth. Typi- 
cally, two such half-power directions exist, ¢; and ¢2, where ¢2 > 1. 
Their difference, A = (¢2 — ¢1), is taken as the beamwidth. 


Ill. BEAMWIDTH RELATIONS 


Returning to the response function (2) and viewing its dependence 
upon the spatial coordinate, we can examine the behavior of H(@) in 
the vicinity of @¢ = ¢’. This corresponds to a value T' = 0, and hence 
to (w7'/2) = 0 in the argument of (2). For small magnitudes of the 
argument we may approximate the denominator term of (2) as 


sin(wT/2) = (wT/2), |wT/2| «K 7/2. (7) 
As a consequence, 


sin{(2N + Mat (8) 


H($) = (2N + 1) | [(2N + 1)(wT/2)] 


The half-power output occurs for a value of the argument of (8) that 
makes the (sin x/x) term equal to 0.707, which requires 


[((2N + 1)(wT/2)] = +1.4 radians. (9) 
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Equation (9) therefore implies the half-power wave-arrival direction 
as 


2.8¢ 
sap = COS cos o’ + (ON + | (10) 


Typically, there will be two half-power directions, ¢; and ¢2, corre- 
sponding respectively to the + and — signs in (10). It is clear from 
(10) that the bracketed term must have a magnitude no greater than 
unity to produce a valid value for ¢-343. As the steering direction is 
made acute to the axis of the array, the value of | cos ¢’ | — 1, which 
indicates a limiting steering direction beyond which (for a given 
frequency and array design) the array will cease to exhibit two direc- 
tions (in the hemisphere) where its output drops to half power. This 
limiting direction, ¢’ = ¢jim, also means that one value of ¢-3a3 must 
be the axial direction (for example, ¢; = 0°, or ¢2 = 180°), which 
implies 


2.8¢ 
fn = cos"? +] 1 — ————_ ], 11 
eee (QN + | on 
The positive sign for the bracketed term corresponds to ¢; = 0° 
with ¢iim in the first quadrant, and the negative sign corresponds to 
2 = 180° with ¢iim in the second quadrant. 
Furthermore, the beamwidth for steering to these limits is simply 


2(2.8c) | 


(2N + l)wd oo 


Aim = cos: — 
Recall that the original approximation requires |(w7')/2| « 7/2, 
which holds for ¢ near ¢’, and for frequencies substantially lower than 
( fupper/2) for any ¢, ¢’. 

The spatial directivity pattern for the one-dimensional array is a 
figure of revolution about the array axis. (In contrast, a two-dimen- 
sional array produces a cigar-shaped beam, confined in two dimen- 
sions.) Equations (10), (11), and (12) describe the conditions for 
finding two half-power beamwidth directions (in the hemisphere 0 < 
@ < 180°) for steering directions no more acute to the axis than jim. 
The complete (figure of revolution) spatial directivity pattern is there- 
fore cone shaped, with the interior of the cone remaining hollow. 
These limits may be particularly appropriate to the case of the array 
mounted on an infinite baffle, such as the wall of a room. 

If the array is unbaffled and freely suspended, steering to angles 
more acute than ¢j;, may be useful. In the limit ¢’ equals 0° or 180°, 
which is the end-fire condition. For all steering directions more acute 
than dim, only one half-power direction exists in the hemisphere, and 
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the interior of the cone-shaped directivity pattern “fills in” until the 
directivity pattern merges into a single main lobe that is axially 
directed. In this case the half-power beamwidth is simply 


A=2¢2, O° <¢<90°, ©’ < dim, 
or 

A = 2(180°-—¢,), 90° <¢@<180°", ’>djim. (18) 
Additionally, the end-fire beamwidth (¢’ = 0° or 180°) is 


} (14) 
IV. REPRESENTATIVE DESIGNS 


To fix the utility of these relations, let us consider a practical set of 
conditions. Suppose the bandwidth of interest is approximately the 
telephone bandwidth. This bandwidth is included within a range 
which, for example, we take as 170 to 3500 Hz. Further, suppose the 
desired range of steering is 60° < ¢’ < 120°, or +30° from broadside. 
For this condition | cos ¢ — cos ¢’ |max = 1.5 and the relation 


fupper = 1/| T'| max = 8500 Hz 


2.8¢ 


(2N + 1)wd 








Aend = 2 feos" 


suggests 
d = 6.5 cm. 


Similarly, the relation 
1M fee ) 
N=|=(7>——-1 
E (; lower nhi 


N = 10. 


These values can be used in eqs. (2) and (10) to examine the spatial 
discrimination and beamwidth under conditions of steering, both 
inside and outside the prescribed range of $’. Figures 4 and 5 illustrate, 
on linear plots, the values of | H(¢)| versus ¢ for ¢’ = 90, 60, 30, and 
0 degrees, and for frequencies two octaves apart, namely, 500 and 2000 
Hz. As mentioned previously, the complete pattern of spatial selectiv- 
ity for the one-dimensional array is a figure of revolution about the 
array axis. 

Consistent with (10), Figs. 4 and 5 show that the acuity of the 
steered beam diminishes with frequency, and with steering directions 
away from the normal to the array axis. In addition, using (11) through 
(14) provides a summary characterization of A versus ¢’, with fre- 


requires 
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STEERING = 90 DEGREES STEERING = 60 DEGREES 





STEERING = 30 DEGREES STEERING = 0 DEGREES 


FREQUENCY = 500 Hz 
(b) 


180 


Fig. 4—Spatial response for the uniform array. The parameters are N = 10, d = 6.5 
cm, and f = 500 Hz for (a) ¢’ = 90°, 60° and (b) ¢’ = 30°, 0°. The polar plots show 
| H()| versus ¢ on a linear amplitude scale. For the one-dimensional array, the complete 
spatial directivity pattern is a figure of revolution about the vertical axis. 
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STEERING = 90 DEGREES STEERING = 60 DEGREES 





STEERING = 30 DEGREES STEERING = 0 DEGREES 


FREQUENCY = 2000 Hz 


(b) 


180 


Fig. 5—Spatial response on a linear scale for the uniform array. The parameters are 
N = 10, d = 6.5 cm and f = 2000 Hz for (a) ¢’ = 90°, 60° and (b) ¢’ = 30°, 0°. 
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HALF-POWER BEAMWIDTH A IN DEGREES 





0 15 30 45 60 75 90 
END-ON BROADSIDE 
STEERING ANGLE @’ IN DEGREES 


Fig. 6—Variation of half-power beamwidth A as a function of steering angle ¢’. 
Frequency is the parameter. The limiting steering angle, jim, is calculated from eq. 
(11). For steering angles more acute to the array axis than ¢iim, only a single half-power 
response direction exists in the right hemisphere. The beam pattern merges to an axially 
directed lobe, as specified in eq. (14), and as indicated by the dashed portions of the 
curves. 


quency the parameter, as shown in Fig. 6. The plot is shown for the 
first quadrant, 0° < ¢’ < 90°. The second quadrant, 90° < ¢’ < 180", 
exhibits mirror-image curves, symmetric about ¢’ = 90°. 

Another aspect of interest is the asymmetry of the beam about the 
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1.50 


ip’ 61)/(G2-6) 
: 


BEAM ASYMMETRY 
Oo 


R= 





0 15 30 45 60 75 90 
END-ON BROADSIDE 
STEERING ANGLE @’ IN DEGREES 


Fig. 7—Variation of beam asymmetry R as a function of steering angle ¢’, with 
frequency as the parameter. 


steering direction $’. This is reflected by the values (¢’ — ¢,) and 
(@2 — o’), and by their ratios, R. This asymmetry, as a function of the 
steering angle and with frequency as a parameter, is plotted in Fig. 7. 


V. CONCLUDING COMMENT 


Automatic, electronic steering of microphone arrays is an attractive 
means for improving sound pickup in large rooms. This improvement 
makes more feasible interactive communication between sizable 
groups of people seated in auditoria or meeting rooms remote from 
one another. Recently emerging microprocessors have the computa- 
tional capability for automatic detection of desired sound sources in 
the room, and for automatic beam forming and beam steering to the 
desired source.’ But spatial acuity is affected by steering, typically 
diminishing as the beam is steered away from the broadside position. 
Therefore, in addition to the considerations of useful bandwidth (and, 
hence, array geometry), the design of a beam-steered system must take 
into account the variations in acuity with steering. The relations 
derived here provide quantitative means for designing uniform, un- 
weighted arrays for specified conditions of room coverage. 
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